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Voorwoord 



Bij liet niaken van een wetenscliappelijk werk is liet steeds belangrijk om er plezier in 
te blijven liebben. Dit was zeker het geval bij deze thesis. Het onderwerp interesseert 
niij zeer sterk, onidat het twee van niijn grote passies conibineert, namehjk wiskunde 
en niuziek. Ik vond het enorni boeiend oni ermee bezig te zijn. 

Een thesis wordt natuurlijk nooit alleen geschreven. Ik wil daarom ook graag mijn 
proniotor Kenneth Sorensen bedanken voor het nalezen en het beantwoorden van 
mijn vragen. Ook mijn vriendin Els is een enorme steun geweest en een hulp bij het 
nalezen. Tenslotte wil ik ook mijn ouders bedanken, die al jaren achter mij staan 
en mij de mogelijkheid hebben gegeven om handelsingenieur te studeren. 
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Inleiding 



Zal IBM de volgende Mozart ontwikkelen? Deze vraag lijkt ver van ons bed te 
staan. Er is echter heel wat aan de gang in de wereld van 'Computer Aided Com- 
posing', of kortweg CAC. Dit relatief nieuwe onderzoeksveld houdt zich bezig met 
het componeren van muziek door computers. 

Wat maakt een muziekstuk beter dan een ander? Dit is een belangrijke vraag die we 
ons moeten stellen wanneer we het gebied van CAC verkennen. Er bestaat zoiets als 
'muziektheorie', deze klassieke wetenschap legt verschihende regels op die bepalen 
wat een goed muziekstuk is. Wanneer we deze regels kunnen kwantificeren en er 
een belangrijkheid aan toekennen, dan krijgen we een soort van 'doelfunctie'. Een 
goed muziekstuk komt overeen met een hogere waarde van zijn doelfunctie. Dit 
doet ons sterk denken aan een optimalisatieprobleem. En inderdaad, er zijn reeds 
verschihende CAC-toepassingen ontwikkeld die gebruik maken van technieken uit 
de wereld van optimalisatie, zoals bijvoorbeeld genetische algoritmen. 

In deze thesis gaan we enkele optimalisatietechnieken, namelijk metaheuristieken, 
meer in detail bekijken. We gaan ze vervolgens koppelen aan CAC. Hoe werkt 
een programma dat muziek kan componeren? Na enkele van deze voorbeelden te 
hebben besproken gaan we zelf een kleine toepassing ontwikkelen, gebaseerd op tabu 
search. We kozen voor tabu search omdat er in het verleden reeds veel projecten 
gebruik maakten van een genetisch algoritme. Van tabu search zijn er weinig of geen 
toepassingen bekend. We hopen hiermee aan te tonen dat een tabu search algoritme 
een nuttig instrument kan zijn bij het componeren van muziek. 

In het eerste hoofdstuk gaan we bespreken hoe we het componeren van muziek 
kunnen beschouwen als een optimalisatieprobleem. We vertellen kort lets meer over 
muziektheorie en wat de 'optimalisatie van muziek' precies inhoudt. Vervolgens 
gaan we dieper in op de geschiedenis van computermuziek. Het lijkt voor de hand 
te liggen dat computers goed overweg kunnen met muziek, maar zoals zal blijken is 
dit niet altijd het geval geweest. 

1 



Inleiding 2 

In hoofdstuk 2 vertellen we lets nieer over de verschillende metalieuristieken. Dit zijn 
combinatorische optimalisatietechnieken die kunnen worden aangewend oni com- 
plexe problenien op te lessen. De volgende nietaheuristieken worden besproken: 
iterated local search, simulated annealing, genetische algoritmen en tabu search. 

Hoofdstuk 3 koppelt de twee vorige hoofdstukken. We bespreken enkele concre- 
te zaken zoals: hoe stellen we een muziekstuk voor in een computerprogramma? 
Welke programmeertalen zijn geschikt om CAC-toepassingen te schrijven? Daaren- 
boven bespreken we enkele voorbeelden van CAC-prograninia's, nanielijk GenJam 
en SIEDP. 

Het tweede deel van deze thesis gaat over de toepassing die we zelf ontworpen hebben 
(TabuMusic). Het ontworpen programma werkt op basis van een tabu search en 
optinialiseert willekeurige honiofone niuziek. De working wordt in detail uiteengezet 
in hoofdstuk 4. In hoofdstuk 5 maken we tenslotte een evaluatie van de resultaten. 

De bijgevoegde cd-rom bevat de broncode van TabuMusic, sanien met een gecompi- 
leerde versie voor Win32 en linux. Het opensource programma iabc is ook bijgevoegd. 
Verder staan er ook enkele geoptimaliseerde muziekfragmenten op. 



DEELI 



Theoretische Achtergrond 
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Muziek als optimaliserings- 
probleem 



/ dream of musical instrum,ents obedient to my thought and which with 
their contribution of a whole new world of unsuspected sounds, will lend 
themselves to the exigencies of my inner rhythm. 

Varese, 1937 zoals aangehaald in (Stephen, 2004) 



1.1 Probleemstelling 

Kunnen computers muziek componeren? Ja, dat kunnen ze zeker, we moe- 
ten de vraag misschien herformuleren als: "Kan een computer met behulp van 
optimalisatie-algoritmen een goed muziekstuk componeren?" 

Eerst en vooral moeten we vaststellen dat een muziekstuk perfect voor te stellen is 
met behulp van een computer. Denk hiervoor maar aan digitale piano's, drumcom- 
puters, sequencers en allerlei andere MIDI-toestellen^. Het gebruik van dit soort 
geluidssynthese is sterk ingeburgerd. 

Hoewel muziek duidelijk is geassocieerd met emoties en artistieke expressie speelt 
logica er een niet te verwaarlozen rol in. Klassieke muziekstukken (en andere) zitten 
zeer logisch en wiskundig ineen en volgen over het algemeen de regels van de har- 
monie (Geeurickx, 1985). Een computer is in staat om met behulp van deze logica 
en regels zelf een muziekstuk te componeren. 



^MIDI staat voor MIDI: Musical Instrument Digital Interface, een gestandaardiseerde inter- 
face gebruikt bij elektronische muziek- keyboards en PC's voor computer controle van muzikale 
instrumenten en apparaten (Disctronics, 2004). In sectie 1.4.3 wordt hierover meer gezegd. 
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Bestaat er een optimaal niuziekstuk? Aan welke criteria moet muziek liiervoor 
voldoen? Wanneer we erin zouden slagen om objectieve criteria in verband met 
harmonie, melodie en ritme vast te leggen, dan kunnen we op basis hiervan een 
doelfunctie opstellen. Op die manier zijn we in staat om een 'score' toe te kennen 
aan muziek, die weergeeft hoe good ze aan de criteria voldoet. Het is dan theoretisch 
mogelijk om een optimaal niuziekstuk te 'berekenen'. Dit is dan een niuziekstuk 
waarbij de doelfunctie zijn maximale waarde bereikt. In de praktijk blijkt het ecliter 
zeer moeilijk om deze criteria ook echt objectief vast te leggen (Kleeven en Philippi, 
2002). 

In de rest van dit hoofdstuk gaan we dieper in op het verband tussen muziek en 
optimalisatie. Verder bespreken we ook de geschiedenis van computermuziek en de 
begrippen die hierbij van belang zijn. 

1.2 Muziek en optimalisatie 

1.2.1 Definitie van muziek 

We konien in het dagelijkse leven bijna elke dag in contact met muziek. Als we over 
straat lopen, in de auto,. . . zelfs op het werk. Maar hoe wordt muziek gedefinieerd? 
Als we de term muziek opzoeken in de Encarta Encyclopedie, dan vinden we volgende 
definitie: 

Music is the organized movement of sounds through a continuum of time 
(Encarta, 1997). 

Muziek is dus een opeenvolging van geluiden doorheen de tijd. Onze ervaring leert 
ons dat niet alle geluiden zoniaar als muziek worden beschouwd. In de volgende sec- 
tie gaan we dieper in op de wiskundige aard van muziek. Waaroni klinken bepaalde 
geluiden good sanien en andere niet? 

1.2.2 Een wiskundige kijk op muziek 

Muziek bestaat uit geluid. Geluiden zijn opgemaakt uit herhalende geluidstrillin- 
gen. De frequentie van deze trillingen wordt gemeten in Hertz (Hz), dit is het aantal 
trillingen per seconde (Land, 2003). In figuur 1.1 is de pianotoon C (129Hz) weer- 
gegeven. De luchtdruk is hier uitgezet tegen de tijd. In de figuur correspondeert de 



1.2 Muziek en optimalisatie 




Figuur 1.1: Grafische weergave van een pianotoon C van 129Hz (Biles, 1994) 

toonhoogte met de frequentie, de luidheid met de amplitude en de klankkleur met 
de specifieke vorm van liet golfpatroon. 

In muziek zijn er slechts twee constante waarden. De eerste is de frequentie van de 
A (la) en bedraagt 440 Hz. De tweede constante waarde is de ratio van de frequen- 
ties tussen een halve toon, deze is v2 = 1.0594630943593. Met beliulp van deze 
twee waarden kunnen de frequenties van alle mogelijke noten worden uitgerekend 
(Blackburn, 2004). Bijvoorbeeld: 

Atj (1 halve toon hoger dan A)= UQHz x v^2 = 466.16 

Pythagoras was de eerste persoon die een verband zag tussen wiskunde en muziek. 
Hij was er sterk van overtuigd dat wiskunde overal was en dat elke waarde als een 
verhouding uit te drukken was. In de tijd van de Grieken was muziek nog niet zo 
ingewikkeld als vandaag en bestonden er maar vijf noten. Pythagoras wees erop dat 
elke noot kon gespeeld worden als een fractie van een snaar (Cox, 2004). 

Stel dat een snaar klinkt als een A. Als we | van de lengte van de snaar nemen, dan 
vinden we de volgende noot (fa C). Verder vinden we: | « D, | fa i?, | ?s F en 
2 ~ yl. Dit systeem wordt voorgesteld in figuur 1.2. Op deze manier kunnen we de 
vijf noten van de Grieken spelen. Men refereert hier soms naar als het 'Pythagoras- 
toonsysteem' (van de Craats en Takens, 2001). 

Hoe zijn we dan aan een systeem met twaalf noten geraakt? Men heeft de ratio's 
van Pythagoras op de andere noten toegepast. Zo vinden we bijvoorbeeld de noot 
B door I de lengte van een E-snaar te nemen. Op deze manier vinden we twaalf 
noten. We kunnen dit grafisch weergeven met een spiraal, zie figuur 1.3. Langs de 
rand is er een schaalverdeling in graden. Deze schaalverdeling laat een logaritmische 
voorstelling toe. 

Een rondgang langs de linkse spiraal (360°) correspondeert met een stijging over 
zeven octaven. Een octaaf meet dan ^r- graden. Een kwint, frequentieverhouding 
2:3, meet dan °^ ^ x ^ « 30.083786 graden. De spaken zijn om de 30 graden 
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Figimr 1.2: De ratio's van Pythagoras 

getekend. De spiraal aan de rechterkant in figuur 1.3 is strakker opgewonden. Een 
rondgang rond de spiraal stelt nu een octaaf voor (in plaats van zeven). Daardoor 
is een kwint nu 210.5865 graden (van de Craats en Takens, 2001). 

Er was ecliter een probleem met dit systeeni. De frequenties van de noten die gevon- 
den werden vertrekkend van een A waren niet helemaal identiek als deze gevonden 
vertrekkend van de Ajj. Dit zorgde voor problenien wanneer twee verschillend ge- 
stemde instrumenten samen speelden (Cox, 2004). Vooral de grote tertsen (F-A, 
C-E en G-B) klinken vals. Neem bijvoorbeeld het interval F-A. In het systeem 
van Pythagoras kunnen we van de F naar de A gaan in vier intervallen (kwinten): 
F-C-G-D-A. Aangezien een kwint een interval is met verhouding g, verhouden de 
frequenties van F en A zich als '^ = ^. Wanneer we de F twee octaven^ hoger 
nemen wordt deze verhouding ^^^ = ^. Wanneer we van dezelfde F vertrekken 
en we nemen de noot die een zuivere grote terts hoger ligt (dit is weer de A), dan 
moeten hun frequenties zich volgens Pythagoras verhouden als 5 = gg- We merken 
op dat er een verschil zit tussen de eerste valse (|y) en de tweede zuivere (|^) grote 
terts (van de Craats en Takens, 2001). Dit kleine verschil in frequentie zorgt ervoor 
dat de geluidstrillingen gaan zweven en daardoor vals klinken (van de Craats en 
Takens, 2001). Figuur 1.4 geeft dit nog eens grafisch weer. 



^Een octaaf hoger komt overeen met een verdubbeling in de frequentie 
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Figuur 1.3: De kwintenspiraal (van de Craats en Takens, 2001) 




Figuur 1.4: De valse grote terts Y-A^w en de zuivere grote terts F-A^^ (van de 
Craats en Takens, 2001) 
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Figuur 1.5: De middentoonspiraal (van de Craats en Takens, 2001) 

De middentoonstemming kwani met een oplossing voor dit probleem. Ze zorgde 
ervoor dat men gebruik maakte van de hierboven vernoemde constante ( v 2) om de 
verhoudingen tussen de noten te bepalen. Op die manier klinken niet enkele noten 
heel vals, maar alle noten een heel klein beetje. We kunnen dit weer in spiraalvorm 
voorstellen, zie figuur 1.5. Zoals op de figuur blijkt is ook de middentoonstemming 
niet perfect. We zien dat Cis (Cjj) en Des (Db) lichtjes van elkaar verschillen. Als 
we naar een pianoklavier kijken, dan valt editor op dat er maar een toets is voor 
Cis en Des. Men zal dus moeten kiezen welk segment van 12 opeenvolgende noten 
men gebruikt uit de middentoonspiraal. In oudere muziek word vaak gekozen voor 
het segment: 

Eb-Bb-F-C-G-D-A-E-B-Ftt-Ctt-Gtj 

De kwint Ab-Eb klinkt dan echter heel erg vals, hij is bijna een kwart toon te klein. 
De klank van de zwevingen die dit interval veroorzaken word vroeger vergeleken 
met het huilen van een wolf, vandaar de naam wolfskwint (van de Craats en Takens, 
2001). 

Deze problemen worden opgelost door de 3 1-tonenstemming van Huygens, zie figuur 
1.6. Hij verbeterde de middentoonstemming door een veel fijnmaziger octaafverde- 
ling toe te passen. Het resultaat is een wiskundig mooie oplossing, waarbij het octaaf 
in 31 delen is verdeeld. In figuur 1.6 zien we dit grafisch voorgesteld. De zuiverheid 
van de geproduceerde muziek zou veel beter zijn als dit systeem zou worden toe- 
gepast. In de praktijk is het echter niet zo gemakkelijk toe te passen, een normale 
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Figuur 1.6: De niiddentoonspiraal en het 31-tonensysteeni van Huygens (van de 
Craats en Takens, 2001) 



piano zou al snel 227 toetsen krijgen. In de hedendaagse muziek wordt daarom nog 
steeds de middentoonstemming gehanteerd (van de Craats en Takens, 2001). 

Op basis van de frequenties van noten zuUen sonimigen beter bij elkaar passen dan 
andere (Bilotta, 2000). Een gevolg hiervan is dat er toonladders^ ontstaan. Be- 
paalde conibinaties van noten zullen 'droevige' muziek produceren, de zogenaamde 
'mineur'-toonladders. In tegenstelling tot de 'diatonische' of 'majeur'-toonladders 
(Schell, 2002). 

Een andere belangrijke consideratie zijn de 'harmonische tonen'. Veel muzikanten 
vergeten dat wanneer we bijvoorbeeld een noot spelen op een viool, we niet alleen 
deze noot horen, maar ook een aantal harmonische tonen. Dit zijn noten met een 
proportionele frequentie aan de originele noot (Cox, 2004). Een voorbeeld maakt 
dit duidelijk: 

Stel dat we een A spelen van 110 Hz. De harmonische tonen die hier samen mee 
klinken zijn voorgesteld in label 1.1. Deze tonen bepalen welke noten een good 
akkoord vormen. Wanneer twee noten tegelijk worden gespeeld en hun harmonische 
tonen versterken elkaar, dan klinkt het akkoord aangenaam. 

We hebben een korte beschrijving gegeven van de onderliggende aard van muziek. De 
regels van de harmonie zijn gebaseerd op hetgeen hierboven is beschreven, zonder dat 
zij daarbij gedetailleerd ingaat op frequenties en ratio's. Uit al deze mathematische 
theorie kunnen we concluderen dat we muziek kunnen beschouwen als geluidsgolven. 



^Een toonladder is een op- of aflopende reeks van tonen in een vaste volgorde (Wikipedia, 2004). 
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Noot 


Frequentie (Hz) 


Grondnoot 


A 


110 


Iste harmonische toon 


A (een octaaf hoger) 


220 


2de harmonische toon 


E 


330 


3de harmonische toon 


A (twee octaven hoger) 


440 


4de harmonische toon 


ctt 


550 


5de harmonische toon 


E 


660 



Tabel 1.1: Harmonische tonen van een A (Cox, 2004) 

en meer bepaald als frequenties. Het lijkt nu duidelijk dat een hedendaagse computer 
geen moeite heeft met het manipuleren van zo'n simpele trillingen. De volgende 
sectie handelt over hoe we muziek kunnen beschouwen als een optimalisatieprobleem. 



1.2.3 Optimalisatie van muziek 

Wat is een optimalisatieprobleem? Als we het componeren van een muziekstuk 
willen bekijken als een optimalisatievraagstuk, dan is het belangrijk dat we eerst de 
eigenschappen hiervan bestuderen. Typisch herkennen we drie basiseigenschappen 
(Dreo et al., 2003): 

1. Er zijn een aantal variabelen 

2. We willen een doelfunctie maximaliseren (of minimaliseren) 

3. Er zijn een aantal beperkingen 

Wanneer we dit vertalen naar het optimaliseren van muziek, dan kunnen we deze 
eigenschappen als volgt invuUen: 

1. Er zijn een aantal variabelen. Welke noten bevinden zich waar in het muziek- 
stuk? Welke duur hebben zij? Hoe lang is het muziekfragment? 

2. Er bestaat een doelfunctie, we kunnen op een bepaalde manier een score toe- 
kennen aan een muziekstuk. Het is onze bedoeling om uiteindelijk het mu- 
ziekstuk met de maximaal (of minimaal, naargelang de berekening) mogelijke 
score te vinden (Towsey et al., 2001). 

3. Kegels uit de muziektheorie zullen een aantal beperkingen opleggen. Zo zuUen 
bepaalde combinaties van noten meer mogen voorkomen dan andere, of zal 
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er een bepaalde structuur in het ritme nioeten zitten. Daarbij komt dat de 
gebruikte noten binnen een bepaald doniein nioeten liggen. Zo kan het domein 
bijvoorbeeld zijn: alle noten die het nienselijk gehoor kan waarnenien of alle 
noten die door een bepaald instrument gespeeld kunnen worden (Towsey et al., 
2001). 

Zoals we reeds hebben opgenierkt is het concept 'optiniale muziek' sterk afhankelijk 
van de manier waarop we een score toekennen. De gebruikte doelfunctie zal meestal 
bestaan uit volgende delen (Kleeven en Philippi, 2002): 

• harmonische samenhang. Hoe good klinken simultane noten samen? 

• melodische samenhang. Hoe good passen opeenvolgende noten bij elkaar? 

• ritmische samenhang. Is het ritme consistent? 

Er is editor geen universele manier ontwikkeld om dit te berekenen, en dat zal 
ook niet gebeuren in de toekomst. In muziek zit steeds een soort van persoonlij- 
ke inbreng. Deze inbreng kunnen we zeer moeilijk met de huidige muziektheorie 
beschrijven (Kleeven en Philippi, 2002). Bepaalde fragmenten zuUen beter klinken 
wanneer ze niet aan de vooropgestelde regels voldoen. Dit kan te wijten zijn aan 
het feit dat de opgestelde regels nog niet voUedig zijn. Het geven van een score aan 
muziek is op dit moment nog steeds zeer subjectief, en brengt veel discussies teweeg 
op het gebied van CAC. 

De berekening van de doelfunctie zal ook sterk verschillen naargelang de muziekstijl 
(Brown, 2002). Dit komt doordat elke muziekstijl zijn eigen typische regels, ritmes 
en melodien heeft. Wij zijn in staat verschillen tussen stijlen te herkennen, maar 
hoe doet een computer dat? Denk aan het volgende voorbeeld: 

Stel dat een computer geprogrammeerd is met drie verschillende doelfuncties: een 
voor klassieke muziek, een voor popmuziek en een voor oosterse muziek. Wanneer 
we nu een klassiek muziekstuk ingeven en we vragen om de score te berekenen met 
behulp van de doelfunctie voor popmuziek zuUen we hoogstwaarschijnlijk een slechte 
score krijgen. Terwijl het stuk misschien een meesterwerk uit de 19de eeuw is, bleek 
het ritme misschien te eentonig, te traag en de melodie te weinig samenhangend 
naar de normen van de popwereld. Met de doelfunctie voor oosterse muziek zouden 
we weer een slechte score krijgen. Wanneer we de doelfunctie voor klassieke muziek 
gebruiken, dan kunnen we pas een waarheidsgetrouwe score kunnen toekennen en 
het stuk evalueren. 
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In het verdere verloop van deze thesis zullen we het vooral hebben over de westerse 
klassieke muziek, tenzij het anders is aangegeven. 

Wat betreft het 'optimaliseren' van muziek moeten we nog aangeven dat het niet 
onze bedoeling is om per se de beste oplossing te vinden. Ondermeer omdat er in 
de vorige secties reeds is aangegeven dat het helemaal niet zeker is dat er zo'n op- 
lossing bestaat. We willen aan de hand van bepaalde beperkingen een redelijk good 
muziekfragment bekomen. In de volgende sectie vertellen we iets meer over een on- 
derzoeksveld dat volop in bloei is, namelijk dat van 'Computer Assisted Composing' 
of CAC. 



1.3 Computer Assisted Composing 

Computer Assisted Composing (CAC) is een nieuw vakgebied dat de laatste jaren 
steeds belangrijker is geworden (Daubresse en Assayag, 2000). Het houdt zich bezig 
met het componeren van muziek met behulp van computers. We kunnen CAC 
opsplitsen in twee fundamentele delen (ThinkQuest, 1996). 

1. Computers gebruiken om het componeren gemakkelijker te maken, als hulp- 
middel. 

2. Computers hun eigen muziek laten componeren. 

In de eerste hetekenis is de computer slechts een hulpmiddel (zij het dan een zeer 
krachtig) in het compositie-proces. Zo kan een componist bijvoorbeeld een com- 
puterprogramma gebruiken dat nagaat of zijn muziekstuk wel speelbaar is. Denk 
hierbij bijvoorbeeld aan een fragment voor gitaar. Het moot mogelijk zijn om de 
akkoorden die worden voorgesteld door de componist fysisch te spelen op het instru- 
ment. Daarvoor zal de computer alle mogelijke vingerposities moeten berekenen en 
controleren. Het zelfde probleem stelt zich ook bij andere instrumenten, zoals bij- 
voorbeeld de piano (Viana en de Morais Junior, 2003). In sectie 3.4 wordt hiervan 
een voorbeeld behandeld, namelijk SIEDP. 

In de tweede hetekenis wordt de computer zelf componist. Het is vooral deze bete- 
kenis die in de literatuur op dit ogenblik het meeste aandacht krijgt. Het principe 
van deze benadering kunnen we meestal in het volgende raamwerk passen (Jacob, 
1995). 

1. Geef criteria op die bepalen wat een good muziekstuk is. 
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2. Bereken een zo optimaal mogelijk niuziekstuk. 

Een niet te onderschatten moeilijkheid ligt bij de keuze van de berekeningswijze. 
Hoe vinden we een goed fragment wanneer er zoveel niogelijke conibinaties zijn. 
We kunnen niet zoniaar alle niogelijke niuziekstukken proberen en zien welke het 
beste is. In dat geval zouden we enorni lang bezig zijn. Er bestaan betere en 
nieer gestructureerde niethodes oni gericht naar een beter muziekfragment te zoeken, 
bijvoorbeeld genetische algoritmen (Jacob, 1995) en cellular automata (Miranda, 
2001). Hierop komen we later nog uitgebreid terug wanneer we het project GenJam 
bespreken in hoofdstuk 3. 

We zuUen ons vooral focussen op de tweede betekenis van CAC: het gebruik van 
de computer om muziek te componeren. De nadruk ligt in deze thesis vooral op 
de algoritmen die gebruikt worden om een goed niuziekstuk te vinden. Het is niet 
onze bedoeling om een cursus muziektheorie samen te vatten, tocli is het interessant 
om kort iets te zeggen over de structuur van muziek en een scliets te geven van de 
geschiedenis van computermuziek. 



1.4 De geschiedenis van computermuziek 

1.4.1 Vroege geschiedenis van klassieke muziek 

Om een algemeen beeld te krijgen van het ontstaan van klassieke muziek, schetsen 
we hier de evolutie van klassieke muziek zonder al te veel in details op te gaan. We 
beginnen in de oudheid en eindigen in het heden. 

In de oudheid begint het verhaal van klassieke muziek. De belangrijkste invloed 
van de oudheid ligt op het gebied van de muziektheorie. Ondermeer de ratio's van 
Pythagoras zitten hier voor iets tussen. Dit is eerder besproken in sectie 1.2.2 
(Wikipedia, 2004). 

Vanaf de vroege Middeleeuwen (tot 1000) is het vooral de kerk die een grote 
invloed uitoefent op de klassieke muziek. De melodieen gezongen in de kerk waren 
eenstemmig en liadden vaak een Aziatische oorsprong (Zychowicz, 2005). 

In de Middeleeuwen (1000-1450) is meerstemmigheid, ook wel polyfonie genoemd, 
een belangrijke vernieuwing. Twee of meer instrumenten spelen hierbij tegelijk. 
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Figuur 1.7: Muziekstuk gegenereerd met het dobbelsteenspel van Mozart (Boskamp, 
1997) 



Geleidelijk aan werd een systeeni van muzieknotatie ontwikkeld waarbij de noot 
wordt voorgesteld als een punt op een balk met lijnen (Zychowicz, 2005). 

Tijdens de Renaissance (1450-1600) werd naast religieuze muziek steeds meer pro- 
fane en instrumentale muziek gemaakt. Er komen strengere regels omtrent de cor- 
recte samenklank, ofwel consonantie van noten. Tenslotte gaat men meer aandacht 
besteden aan de relatie tekst-muziek (Wikipedia, 2004). 

De Barok (1600-1750) zorgde voor een verandering van muziekstijl op vijf jaar 
tijd^. De meeste moderne muziekinstrumenten werden ontwikkeld in de Barok: de 
strijkinstrumenten en blaasinstrumenten. De Barok eindigt met de dood van de 
componist Johann Sebastian Bach (1685-1750) (SonyMusicEntertainement, 2001). 

Het Classicisme (1750-1810) geeft liaar naam aan de klassieke muziek. Deze korte 
periode wordt gekenmerkt door componisten zoals Wolfgang Amadous Mozart en 
Joseph Haydn. Nieuwe vormen van muziek zoals de sonatevorm en de symfonie 
ontstaan in het Classicisme (Wikipedia, 2004). 

De componist Mozart heeft in deze periode reeds een bijdrage geleverd aan de wereld 
van het CAC. Mozart is namelijk de uitvinder van het 'musikalisches wlirfelspiel', 
ofwel 'muzikaal dobbelsteenspel'. Dit spel laat toe om willekeurige muziekstukken 
te genereren met behulp van twee dobbelstenen en 12 korte muziekfragmenten (met 



^Dit komt ondermeer omdat de monodie, met het systeeni van 'basso continuo ' opkomen, sanien 
met de harmonie. 
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als lengte 1 maat) (Burns, 2005). Zijn idee was dat de spelers in staat waren een 
groot aantal nielodien te genereren door willekeurig te selecteren uit een aantal 
keuzemogelijkheden (Peterson, 2001). Voor het componeren van een niuziekstuk 
van 16 maten gaat het spel als volgt (MathsoftEngineering en Education, 2004): 

Doe voor elke maat m van 1 tot 15 (exclusief m = 8) het volgende: 

1. Gooi de dobbelstenen zodat je een getal i krijgt dat elf mogelijke waarden kan 
aannemen (van 2 tot en met 12). 

2. Kies als opvulling van maat jtt, het fragment dat bij i hoort 

Deze eenvoudige regel geldt voor elke maat, behalve voor maat 8 en 16. Hiervoor 
zijn enkel fragment 2 of 3 mogelijk. De beslissing van welk fragment we nemen zal 
afhangen van het feit of het gegooide getal i even of oneven is. De reden voor deze 
laatste regel ligt in de muziek die we genereren. Typisch klinkt deze vrij sonate- 
achtig en zorgen de 8ste en 16de maat voor een rustpunt (zie figuur 1.7). Deze 
methode op zich is redelijk simpel en produceert toch goede muziek. Een niet te 
verwaarlozen voorwaarde is echter het hebben van de 12 korte fragmenten waarop 
alles is gebaseerd. Mozart heeft er een heleboel achtergelaten en wanneer we het 
spel hierop toepassen krijgen we een verrassend muziekje met een duidelijke 'Mozart- 
tint'. Een mogelijk resultaat van het dobbelsteenspel is weergegeven in figuur 1.7. 

Een volgende periode is de Romantiek (1810-1910). Deze muziek maakt ongeveer 
95% uit van wat op dit moment in concerten te horen is. De Verlichtingsgedach- 
te stimuleerde diverse muzikale ontwikkelingen. Dit leidde tot complexe harmoni- 
sche ontwikkelingen, steeds verbeterende muziekinstrumenten, grotere muziekstuk- 
ken, enzoverder (Zychowicz, 2005). 

Vanaf de 20ste eeuw bestaan er verschillende stromingen gelijktijdig: het Moder- 
nisme, Neo-stijlen zoals het Neo-Calassicisme en Neo-Barok, Avant-Gardise, Seria- 
lisme. Minimal music. New Complexity, Nieuwe Eenvoud en Post Modernisme. We 
beperken ons tot een opsomming van deze stijlen, een voUedige bespreking zou ons 
te ver leiden (Blood, 2004). 

1.4.2 Nieuwe instrumenten en machines 

Muziekinstrumenten bestaan al heel lang. Ze worden in de loop van de geschiede- 
nis voortdurend bijgewerkt en verfijnd. In het begin van de 20ste eeuw beginnen 
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mensen echter een ander soort instrumenten te ontwikkelen, de elektronische mu- 
ziekinstrumenten. In 1906 registreert Thaddeus Caliil in de Verenigde Staten een 
patent op een 'elektrisch gebaseerd geluidsgeneratie systeem', beter gekend als de 
'dynamofoon' of het 'telharnionicuni' (tadream.net, 2005), zie figuur 1.8. Hoewel de 




Figuur 1.8: De dynamofoon (Sonhors, 2005) 

afnietingen van het apparaat indrukwekkend waren (200 ton gewicht en 1.8 meter 
hoog), kreeg de uitvinding niet veel aandaclit. Deze machine was het eerste additie- 
ve geluidssyntheseapparaat en kon geluid produceren van verschillende frequenties 
en intensiteit. In het begin van de 20ste eeuw werden meerdere van deze synthese- 
apparaten gebouwd door onder meer Theremin, Milhaud, Mager en Martenot. Deze 
'instrumenten' waren allemaal in staat om geluiden te creeren die vroeger ondenk- 
baar waren. Helaas waren ze geen groot succes. Een eerste reden is dat ze allemaal 
een beperkte capaciteit hadden en vrij primitief gebouwd waren. Verder waren er 
nog geen goede methodes om de geproduceerde geluiden op te slaan en te bewerken. 
Ten slotte wordt vaak beweerd dat de muziekwereld nog niet klaar was voor instru- 
menten die een radicale verandering teweeg brachten. Vele componisten bleven de 
klassieke Westerse instrumenten trouw (Stephen, 2004). 

In 1947 was er een doorbraak op het gebied van CAC. Lejaren Hiller en Leonard 
Isaacson experimenteerden met de Illiac I'-mainframe, een Univac computer (figuur 
1.9), aan de universiteit van Illinois (Burns, 2005). Met hun onderzoek wilden ze 
aantonen dat (Pearce et al., 2002): 

• muzikale technieken kunnen worden verwerkt door computerprogrammering. 

• computers kunnen worden gebruikt door componisten om hun compositietech- 
nieken uit te breiden. 

• computers op een ongebruikelijke manier kunnen worden ingezet om een radi- 
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Figuur 1.9: De lUiac-mainframe (MediaArtNet, 2004) 



caal ander soort niuziek te bekomen. 

Hun experiment bestond uit vier delen, met toenemende complexiteit (Kleeven en 
Philippi, 2002): 

• experiment 1: generatie van muziek met een beperkt aantal regels. 

• experiment 2: generatie van muziek met waarbij telkens meer regels werden 
toegevoegd. 

• experiment 3: implementatie van ritme, dynamiek en tempo. 

• experiment 4: generatie van muziek op basis van Markov ketens. Aan de 
hand van alle voorgaande noten word de kans berekend dat de volgende noot 
bijvoorbeeld een C zou zijn. De noot met de lioogste kans word gekozen. 

Ze waren de eersten die erin slaagden om een computer een muziekstuk voUedig 
zelf te laten componeren. Het resultaat is de bekende 'lUiac Suite', ook bekend 
als 'String Quartet No. 4', uitgebracht in 1956 (ThinkQuest, 2004). We moeten 
opmerken dat hoewel de muziek word geschreven door de computer, het stuk word 
gespeeld door echte instrumenten. Met andere woorden, de Illiac I was niet in staat 
om de muziek te genereren (MediaArtNet, 2004). 

Op het vlak van elektronische muziek beginnen er echter serieuze veranderingen te 
komen vanaf 1948. Pierre Schaeffer was de eerste die er begon te experimenteren 
met elektronische muziek. Hij nam geluiden op uit het alledaagse leven, bewerkte 
en mixte ze. Zijn 'Musique Concrete' was het eerste muziekstuk dat word opgevoerd 
zonder dat er mensen aan te pas kwamen. Schaeffer kreeg vele volgelingen. Op deze 
manier begon het elektronische muziektijdperk echt (Stephen, 2004). 
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1.4.3 Computers en muziek 

Een van de eerste eclite computerprogramma's voor muziek werd geschreven in de 
jaren '60 door Max Mathews (Stephen, 2004). Matthews werkte bij Bell Labs, een 
onderzoekscentruni dat door AT&T werd gerund. Hij niaakte er voorwerpen voor 
het telefoonbedrijf, zijn vrije tijd gebruikte hij oni muzieksoftware te ontwerpen. De 
software die hij schreef heette 'MUSIC. In het begin was het programma in staat om 
enkele tonen te produceren. De volgende versie, MUSIC II, kon vier tonen tegelijk 
produceren. Op een van de snelste computers van die tijd koste het echter nog 
ongeveer een uur computertijd om twee minuten muziek te genereren. De volgende 
versies, MUSIC IV en V waren efficienter en konden nog meer geluiden produceren 
(tadream.net, 2005). Steeds meer componisten begonnen zich te vertonen in het 
labo van Matthews. Uiteindelijk werden er zelf enkele componisten als assistent 
aangenomen en werd het project een volwaardig onderzoeksproject (UCSC, 2004). 

In de jaren '70 werden er een heleboel onderzoeksgroepen opgericht die computermu- 
ziek bestudeerden. De bekendste waren die van M.I.T., De Universiteit van Illinois 
in Champaigne-Urbana, de Universiteit van California in San Diego en het 'Institut 
de Recherche et Coordination Acoustic/Music' in Parijs (IRCAM) (UCSC, 2004). 
Deze centra werkten allemaal met een enorme mainframe. Een muziekstuk compo- 
neren was toen een langdurig proces. Na enkele dagen tekst ingeven, compileerde de 
computer de geluiden op een cassette of harde schijf. Vervolgens moest dit nog eens 
door een digitaal-analoogconverteerder worden overgezet op een cassette. Dan pas 
kon men de muziek horen. In deze tijd werden er dan ook meer papers geschreven 
dan muziek gecomponeerd (UCSC, 2004). 

In 1968 ontwikkelt Max Mathews een van de eerste hyhride systemen. Dit wil zeggen 
dat hij een koppeling maakte tussen een computer en een synthesizer. Terwijl een 
muzikant op een synthesizer speelde nam het GROOVE programma van Mathews 
alles op. Daarna kon deze muziek bewerkt worden door de computer en terug 
afgespeeld worden door de synthesizer. Andere onderzoekers hebben enkele variaties 
hierop gemaakt (CSounds, 2004). 

Het grote succes van hybride systemen komt samen met de komst van de eerste 
thuiscomputers. De Apple II maakte het bijvoorbeeld bijzonder gemakkelijk om een 
synthesizer aan te sluiten. Een breed publiek kwam nu in contact met de wereld 
van computercompositie. Ook de komst van de Commodore 64, in 1982 was een 
doorbraak. Dit type computer bevatte een vierstemmige synthesizer op een chip 



1.4 De geschiedenis van computermuziek 20 

(UCSC, 2004). 

Tegen 1983 kwam de MIDI-standaard (Musical Interface Digital Instrument) er. 
Deze standaard maakt liet niogelijk oni gelijk welke computer met eender welke syn- 
thesizer te verbinden. De verkoop van MIDI-synthesizers kende een enorme 'boom' 
in de jaren '80 (Glatt, 2004). Het MIDI-systeem is een communicatieprotocol voor 
elektronische muziek dat zowel hardware als software combineert. Een belangrijk 
concept hierbij is een MIDI-event. Een event stelt de ontvangst (of verzending) van 
een aantal bytes voor. Twee veel voorkomende events zijn (Wiggins et al., 1993a): 

• Note On 

• Note Off 

Elke event bestaat steeds uit drie delen (Wiggins et al., 1993a): 

• een identificeerder, of ID 

• een toonhoogte (uitgedrukt in een integer-waarde) 

• het tijdstip waarop de noot wordt gespeeld 

We kunnen al deze informatiestromen bundelen in een bestand, de zogenaamde 
MIDI-file. Hierin zit alle informatie vervat om een fragment te spelen op een syn- 
thesizer of MIDI-instrument. Aangezien we een MIDI-bestand niet in ASCII-vorm 
(tekstvorm) kunnen lezen, hebben we gebruik gemaakt van het conversieprogram- 
ma 'MIDI File DisAssembler'. Een tool die gratis kan afgeladen worden op het net 
(Glatt, 1998). De belangrijkste functie van het programma is dat het een MIDI- 
bestand kan omzetten in een leesbaar ASCII-bestand. De tekstversie van een wille- 
keurig gekozen muziekstuk is gegeven in bijlage A. 

In de tekstuele vorm van het MIDI-bestand zien we de drie eigenschappen van een 
event duidelijk terugkomen: 

• ID: de nummer voor een event 

• toonhoogte of 'pitch': E 4 voor het eerste event -^ een E (mi) in het vierde 
octaaf 

• de snelheid van indrukken: die blijft hier onveranderlijk en bedraagt nul. 

In het begin van het bestand wordt additionele informatie gegeven over het muziek- 
stuk. Bijvoorbeeld de toonaard, het tempo, de maatsoort, . . . 
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In de onderzoekscentra van computermuziek focuste men zich op het ontwikkelen 
van software. De meeste MIDI-instrumenten waren eigenlijk gewoon heel simpele 
computers die een bepaald programma draaiden. Op M.I.T. word een versie van het 
MUSIC programma gemaakt die kon worden gecompileerd op bijna alle platformen. 
De mainframes werden vervangen door desktops en niet lang daarna word de taal 
C Sound ontwikkeld. Dit is een zeer krachtige taal die nu vaak wordt gebruikt voor 
het ontwikkelen van compositie-programma's (Stephen, 2004). Ondertussen zijn er 
al meer talen en toepassingen ontwikkeld (zie hoofdstuk 3. De wereld van CAC 
krijgt steeds meer mogelijkheden . . . 

In het volgende hoofdstuk gaan we dieper in op een bepaalde soort optimalisatie- 
technieken, namelijk metaheuristieken. 



2 



Metaheuristieken 



// you can't solve a problem, then there is an easier problem you can 
solve: find it. 



(Polya, 1957) 



2.1 Context 

Metaheuristiek is een vakgebied in opbloei. Dit is onder meer te wijten aan de steeds 
grotere rekenkracht van computers. Bepaalde metaheuristieken die vroeger niet te 
implementeren waren, kunnen nu gemakkelijk worden verwerkt door computers. Een 
precieze definitie en omschrijving van deze methodes wordt gegeven in sectie 2.3. 
Metaheuristieken zijn optimalisatietechnieken die vaak gebaseerd zijn op processen 
uit andere vakgebieden die niets te maken hebben met optimahsatie, zoals biologie 
en fysica. Dit zal bhjken wanneer we de verschihende metaheuristieken verder gaan 
bespreken (Pirlot, 1996). 

Voor welke problemen zijn metaheuristieken de geschikte oplossingsmethode? Een 
voorstelling van de verschillende optimalisatieproblemen en de manier waarop ze 
worden opgelost wordt gegeven in figuur 2.1. 

In figuur 2.1 vertrekt men van de basiseigenschappen van optimalisatieproblemen: 
ze minimaliseren steeds een doelfunctie, hebben een aantal (doel-)variabelen die 
geidentificeerd moeten worden en hebben een aantal beperkingen. Deze drie eigen- 
schappen karakteriseren het voUedige probleem. We kunnen vervolgens een grote 
opdeling maken tussen discrete en continue problemen. Discrete problemen zijn 
problemen waarbij de uitkomstenverzameling 'telbaar' is, we kunnen de uitkomsten 
als het ware tellen (Willeme en van Rompay, 2001-2002). Deze verzameling kan een 
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Figuur 2.1: Algemene classificatie van optimalisatiemetliodes (Dreo et al., 2003) 



eindig of oneindig aantal elementen bevatten. In het continue geval is de uitkom- 
stenverzanieling onaftelbaar en bevat zij steeds een oneindig aantal elementen (van 
Nuffelen, 1998). 

Bij de discrete problemen interesseren we ons vooral voor de comhinatorische pro- 
blemen. Dit zijn optinialiseringsproblenien die de beste conibinatie van een aantal 
variabelen of variabele-waarden proberen te vinden (Wikipedia, 2004). Een voor- 
beeld van een discreet probleem is het bekende Travelling Salesman Problem (TSP). 
Een liandelsreiziger nioet n steden exact een keer bezoeken. In zijn eenvoudigste 
vorm bestaat het TSP probleem uit het minimaliseren van de totale reisduur, waarbij 
de afstanden tussen de steden gekend zijn. Dit probleem kan nog worden uitgebreid 
door bijvoorbeeld een beperking op de tijd waarop hij een bepaalde klant moot 
bezoeken, het aantal 'folders' dat hij kan meenemen, enz. We komen later nog uit- 
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gebreid op dit voorbeeld terug. In bepaalde gevallen kan dit soort probleeni een 
berekening vergen waarvan de grootte exponentieel toeneenit met het aantal steden. 
De volgende illustratie maakt ons het belang van exponentiele groei duidelijk: 

Een oude legende vertelt het verhaal van de uitvinder van het schaakspel 
en zijn opdrachtgever, de koning van Persie. Naar verluidt was de koning zeer 
tevreden met dit nieuwe spel en daarom wou hij de uitvinder een weelderige 
compensatie geven. De man zou genoeg eten krijgen om zijn familie een jaar 
te voeden. De uitvinder, een wiskundige begaafd man, bedacht een slimme 
manier om meer te krijgen. 

Hij kwam tot een akkoord met de koning. Op de eerste dag zou de 
koning een rijstkorrel op het eerste vakje van het schaakbord ieggen voor de 
uitvinder. De tweede dag zou hij dubbel zoveel rijstkorrels op het tweede 
vakje ieggen. De derde dag dubbel zoveel op het derde vakje, enzovoorts 
tot het 64ste vakje (Bick, 2000). 

De koning vond dit een fantastisch idee en hij Net zijn bedienden de 
eerste dag een korrel rijst klaarleggen voor de uitvinder. Zoals afgesproken 
kwam de man eike dag de klaargelegde rijstkorrels halen om zijn gezin te 
voeden. De eerste dagen noemde iedereen hem een idioot en vreesde zijn 
gezin zelfs voor ondervoeding. 

Maar het tij keerde. Het aantal rijstkorrels die hij mocht komen afhalen 
verdubbelde eIke dag en tegen de negende dag waren het er al 256 (2'='), op 
de tiende 512, op de elfde 1024, 2048 op de twaalfde dag . . . Tegen de 64ste 
dag was de koning meer dan 10 000 000 000 000 000 000 of 10^^ rijstkorrels 
verschuldigd aan de uitvinder. Om zo'n hoeveelheid rijst de groeien is een 
oppervlakte nodig die groter is dan de aarde, inclusief de oceanen. 

De koning kon deze vraag niet voldoen en eindigde als schuldenaar, terwiji 
de uitvinder van het schaakspel en zijn familie hun intrek namen in het paleis 
(Dittii, 2003). 

Zo'n exponentiele groei van de rekentijd willen we zeker vermijden. Om dit te berei- 
ken zijn er benaderende methodes zoals specifieke heuristieken en metaheuristieken. 
Een belangrijke eigenschap van deze technieken is dat ze op een relatief korte tijd 
tot een redelijk goede oplossing kunnen komen, in plaats van jaren te rekenen om 
de allerbeste oplossing te vinden. 

In het continue geval maakt men een onderscheid tussen lineaire en niet-lineaire 
problemen. Lineaire problemen kunnen worden opgelost met behulp van lineair pro- 
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granimeren. Niet lineaire problenien zijn moeilijk op te lossen. We kunnen gebruik 
maken van lokale metliodes, of ze discretiseren en met beliulp van metalieuristieken 
oplossen (Dreo et al., 2003). 

Samenvattend kunnen we stellen dat metaheuristieken zullen worden gebruikt in 
twee gevallen: 

• discreet: wanneer de onivang van het probleem exponentieel toeneenit. We 
kunnen theoretiscli gezien wel tot de beste oplossing komen, niaar dit zal 
enorm veel rekentijd in beslag nenien. 

• continu: wanneer er geen algoritnie gekend is dat naar een globaal optimum 
leidt na een beperkt aantal bewerkingen. 

In de literatuur rond optimalisatie en artificiele intelligentie bestaan er een aantal 
klassieke problemen. Zo is er bijvoorbeeld het al eerder besproken travelling sa- 
lesman problem, het knapsack-probleem, het vehicle routing probleem,. . .We gaan 
deze problemen kort schetsen, zodat we in de volgende hoofdstukken de draad terug 
kunnen oppikken en concrete oplossingsmethodes voorstellen voor een aantal van 
deze complexe problemen. We merken op dat ze slechts metaforen zijn. We kun- 
nen andere problemen met behulp van deze klassieke problemen interpreteren en zo 
gebruik maken van de oplossingsmethoden die worden voorgesteld. 

Travelling salesman problem 

De term 'Travelling Salesman' is voor het eerst gebruikt in 1932 in het Duitse 
book: "The travelling salesman, how and what he should do to get commissions 
and be successful in his business". Dit book is geschreven door een veteraan- 
handelsreiziger (Michalewicz, 1992). In 1948 heeft RAND Corporation dit probleem 
opnieuw geintroduceerd. Nadien is het een zeer bekend en populair probleem ge- 
worden, onder meer door de opkomst van technieken zoals lineair programmeren 
(Michalewicz, 1992). 

Het travelling salesman problem (TSP) handelt, zoals reeds blijkt uit de naam, over 
een handelsreiziger. Deze moot ?2 steden exact een maal bezoeken en dan terugkeren 
naar zijn vertrekpunt. Hoe moot hij zijn reis plannen, gegeven de kost om tussen 
alle steden te reizen (Luger en Stubblefield, 2003)? Dit interesseert ons omdat het 
aantal mogelijke oplossingen exponentieel toeneenit. Het zal dus zeer lang duren om 
de allerbeste volgorde van steden te vinden (Corne et al., 1999). 
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Knapsack problem 

Het 'knapsack problem' of rugzakprobleeni is een ander voorbeeld van een klassiek 
optinialisatieprobleem. We hebben een verzanieling van n items, die allemaal een 
bepaald gewicht en een winst hebben. Het probleem stelt zich wanneer we een 
rugzak met een capaciteit C optimaal willen vullen. Onze doelstelling is het vinden 
van een deelverzameling die optimaal is, waarbij de winst maximaal is. Dit terwijl 
het totale gewicht de capaciteit C niet overschrijdt, zodat de oplossing geldig blijft. 
Een geldige oplossing wordt ook 'feasible' genoemd (Sorensen). Een voorbeeld van 
gegevens voor een knapsack probleem wordt gegeven in label 2.1 (zie verder). Dit 
wordt besproken in sectie 2.3.3. 

Vehicle routing problem 

Bij het 'vehicle routing problem', of simpelweg VRP, hebben we een vloot van voer- 
tuigen gestationeerd in een of meerdere depots. Er zijn ook een aantal geografisch 
verspreide steden, met elk een gekende vraag, die moeten worden bezocht. De op- 
gave is om de verzameling van routes, elk vertrekkend en eindigend in een depot, te 
optimaliseren. Dit wil zeggen dat de kost minimaal moot zijn. In figuur 2.2 wordt 
een voorbeeld gegeven van een mogelijk VRP probleem met een depot (Mastrolilli, 
2005). 

Dit soort problemen komen veel voor in de praktijk. Bijvoorbeeld op het gebied van 
logistiek, bij transportfirma's en in de distributiesector. Zo heeft de firma 'IJsboerke' 
enkele jaren geleden een systeem geimplementeerd dat elke morgen de routes van 
de koelwagens berekende. We herkennen hier duidelijk een VRP systeem (hetgeen 
zij 'Home Vending Support Application noemen'): er zijn een aantal depots waar 
koelwagens kunnen vertrekken, er zijn een aantal klanten met een gekende vraag (of 
toch een voorspelbare vraag) en vertrekkend van deze gegevens wil men de kortste 
route bepalen (Selligent, 2005). 

In de praktijk kunnen er natuurlijk nog extra beperkingen bijkomen. Enkele veel- 
voorkomende worden hier opgesomd (VRP Web, 2005): 

• een beperking op de capaciteit van een voertuig. Bijvoorbeeld: een vracht- 
wagen kan maar maximum 10 ton vervoeren. Dit probleem is gekend als 
'Capacitated VRP' of CVRP (Gendreau et al., 2001). 
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Figuur 2.2: CVRP met een depot, op basis van (Corne et al., 1999) 

• elke klant iiioet bediend worden binnen een bepaald tijdsvenster, ook wel be- 
kend als 'VRP with time windows' of VRPTW (Corne et al., 1999). 

• er zijn meerdere depots, 'Multiple Depot VRP' of MDVRP. 

• klanten kunnen goederen terug meegeven naar het depot, 'VRP with Pick-Up 
and Delivering of VRPPD. 

• klanten kunnen bediend worden door meerdere voertuigen, 'Split Delivery 
VRP' of SDVRP. 



2.2 Heuristiek 



Heuristiek is de kunst van het ontdekken en uitvinden. Het woord heuristiek stamt af 
van het Griekse woord 'heuriskein' (evpiaK.en'), wat 'vinden' betekent (Wikipedia, 
2004). Denk maar aan de kreet die Archimedes riep toen hij uit zijn befaamde 
badkamer kwam, 'Eureka', wat zoveel betekent als 'Ik heb het gevonden'. 
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Figuur 2.3: George Polya (CIS, 2004) 

Een van de personen die meegeholpen heeft oni heuristiek populair te niaken in de 
20ste eeuw is George Polya. In 1957 sclireef liij een book (Polya, 1957), waarin liij 
de verschillende probleemoplossende methodes beschrijft. Hij onischrijft heuristiek 
daarin als: "the study of the methods and rules of discovery and invention" . Als kind 
was hij gefascineerd door de manier waarop wiskundigen bewijzen ontwikkelden. Hij 
leerde op school de bewijzen, maar hij leerde niet hoe deze tot stand zijn gekomen. 
In zijn book stelt hij volgende stappen voor om een wiskundig probleeni op te lessen 
(Wikipedia, 2004): 

1. Eerst nioet je het probleeni begrijpen 

2. Stel vervolgens een plan op 

3. Voer het plan uit 

4. Kijk terug naar je work, wat kan je eraan verbeteren? 

Heuristische algoritnien gebruiken probleemkennis om de oplossingsruimte te door- 
zoeken en zo de explosie hiervan te beperken. Ze proberen het zoeken te richten 
naar de meest belovende regie. Het nadeel van deze heuristische methodes is dat ze 
vaak degenereren naar een exhaustieve zoekfunctie (Silver, 2002). 

Vanaf het begin van de jaren '70 word het duidelijk dat het praktisch onmogelijk is 
om een optimale oplossing te vinden voor problemen zoals VRP en TSP. De rol van 
heuristieken bij het oplossen van combinatorische problemen uit de praktijk word 
steeds belangrijker (Gendreau, 2002). 
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We kunnen dit verduidelijken met het eerder genoenide probleeni van de Travelling 
Salesman. In het uitgewerkte voorbeeld nioet een handelsreiziger negen (A tot I) 
steden bezoeken, vertrekkende en eindigend in stad A. In nieer realistischer proble- 
men zal het aantal steden echter veel hoger zijn. In figuur 2.4 wordt de volgende 
heuristische strategie voorgesteld: 

• Vertrek vanuit stad A naar een wihekeurige stad, bijvoorbeeld stad B 

• Voor elke stad B tot H: 

— Kijk met je aangezicht naar de stad waarvan je komt 

— Begin nu naar rechts te draaien 

— Ga naar de eerste stad die je tegenkomt 




Figuur 2.4: TSP opgelost met een heuristische strategie 

We kunnen uit de grafiek concluderen dat deze snel gevonden oplossing niet de 
optimale zal zijn. Zoals eerder al aangegeven had het veel langer geduurd om de 
allerbeste oplossing voor dit probleem te vinden. Er bestaan echter methodes die 
tot een kwalitatief betere oplossing leiden dan klassieke heuristieken, terwijl ze veel 
minder rekentijd nodig hebben dan exacte methodes, namelijk metaheuristieken 
(Silver, 2002). 
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2.3 Metaheuristiek 

2.3.1 Definitie 

Zoals de naam al laat vernioeden houden metaheuristieken (MH) verband met de 
klassieke lieuristieken. Het woord metaheuristiek bestaat uit twee delen: 'meta' en 
'heuristiek'. De term meta betekent 'verder', 'op een hoger niveau'. Metaheuristie- 
ken zijn geen heuristieken, maar raamwerken. Het zijn richtlijnen om een heuristiek 
te maken. Omdat ze specifiek op een bepaald probleem worden toegepast, zal het 
resulterende algoritme verschillen van probleem tot probleem. Metaheuristieken zijn 
toepasbaar op vrijwel alle combinatorische optimalisatieproblemen (Sorensen). 

We vinden volgende definities terug in de literatuur: 

A metaheuristic is a set of concepts that can be used to define heuris- 
tic methods that can be applied to a wide set of different problems. In 
other words, a metaheuristic can be seen as a general algorithmic fra- 
mework which can be applied to different optimization problems with 
relatively few modifications to make them adapted to a specific problem 
(MetaheuristicsNetwork, 2004) . 

A metaheuristic is formally defined as an iterative generation process 
which guides a subordinate heuristic by combining intelligently different 
concepts for exploring and exploiting the search space, learning strategies 
are used to structure information in order to find efficiently near-optimal 
solutions (Osman en Laporte, 1996), zoals aangehaald in (Blum en Roll, 
2003). 

In tegenstelling tot een exacte oplossingsmethode zoekt een metaheuristiek niet de 
allerbeste oplossing, maar probeert op een snelle manier tot een goede oplossing te 
komen. Deze unieke balans tussen kwaliteit en snelheid zorgt ervoor dat metaheu- 
ristieken steeds meer worden gebruikt voor het oplossen van diverse optimalisatie- 
problemen. Figuur 2.5 geeft de verschillen weer tussen exacte oplossingsmethodes, 
heuristieken en metaheuristieken. 
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Figuur 2.5: Vergelijking van verschillende optimalisatieteclinieken 

2.3.2 Globale indeling 

Er zijn een aantal verschillende metalieuristieken. We kunnen ze in drie grote cate- 
gorien onderverdelen (Blum en Roll, 2003), zoals ook te zien is in figuur 2.1: 

1. local search gebaseerde metalieuristieken 

2. populatiegebaseerde metaheuristieken 

3. hybride metaheuristieken 

In de hierop volgende secties worden deze drie groepen meer in detail besproken. 



Local Search gebaseerde MH 

Bij local search MH gaan we van de ene naar de andere oplossing door kleine wijzi- 
gingen aan te brengen (de zogenaamde moves). Enkele voorbeelden van moves bij 
het handelsreizigersprobleem zijn: 

• een stad op een andere positie zetten (zodat ze bijvoorbeeld later of vroeger 
wordt bezocht) 

• de volgorde van twee steden omdraaien 

• de volgorde van drie steden wisselen 
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De neighbourhood N{xi) is de verzameling van alle oplossingen die we vanuit een 
huidige oplossing (xi) kunnen bereiken door middel van een move (Focacci et al., 
2001). De vraag is natuurlijk welke oplossing uit de neighbourhood we als volgende 
oplossing gaan aanvaarden, X2 G N{xi). Deze beslissing wordt bepaald door de 
strategie die we volgen. Twee voorbeelden van mogelijke strategieen zijn: 

• steepest descent: hierbij aanvaarden we de beste oplossing in de neighbour- 
hood. 

• mildest descent: we genereren de neighbourhood maar stukje voor stukje. 
Vanaf het moment dat we een betere oplossing tegenkomen, dan aanvaarden 
we die (Silver, 2002). 

Het concept van local search wordt nog verder geillustreerd in sectie 2.4.4 over tabu 
search. 

Een vaak voorkomend probleem is dat een strategie blijft vastzitten in een lokaal 
optimum. Een lokaal optimum is een oplossing waarbij onze strategie geen betere 
oplossing meer kan vinden, terwijl er nog wel een betere bestaat. Op de linkerhelft 
van figuur 2.6 is dit grafisch weergegeven. Vanuit punt x zou het kunnen dat de 
strategie niet verder kijkt dan de stijging rond x en dus niet weet dat er nog een 
lager punt is (Blum en Roli, 2003). De kracht van metaheuristieken is dat zij kunnen 
ontsnappen uit zo'n lokaal optimum en daardoor meer kans maken om het globale 
optimum, de allerbeste oplossing, te bereiken, zoals is voorgesteld in figuur 2.6. Het 
maken van een move wordt herhaald totdat een bepaald stop-criterium bereikt is. 



Figuur 2.6: Lokaal optimum versus globaal optimum 
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Enkele bekende local search nietalieuristieken zijn tabu search, simulated annealing, 
iterated local search. Deze worden uitgebreid besproken verder in dit hoofdstuk 
(secties 2.4.1, 2.4.2 en 2.4.4). 

Populatiegebaseerde MH 

Populatiegebaseerde MH hebben hun wortels in de genetische biologie en nemen 
dan ook een deel van deze nonienclatuur over. Ze behandelen steeds een aantal 
oplossingen sanien, deze vormen een populatie. Oni tot een volgende oplossing te 
komen, combineren ze bestaande oplossingen, waarbij hopelijk de beste eigenschap- 
pen van de 'ouders' behouden blijven. Het combineren gebeurt tot een bepaald 
stop-criterium bereikt is. Kenmerkend voor dit soort metaheuristieken is dat zij op 
een vrij willekeurige manier de oplossingsruimte verkennen, in tegenstelling tot local 
search strategieen die meer gericht verkennen (Dreo et al., 2003). 

Bekende voorbeelden van populatiegebaseerde metaheuristieken zijn evolutionaire 
algoritmen, ant colony optimization,. . . Evolutionaire algoritmen worden verder be- 
sproken in sectie 2.4.3. 

Hybride MH 

Als derde categorie kunnen we de hybride metaheuristieken beschouwen. De laatste 
tijd is er een toenemende interesse in de literatuur vast te stellen voor hybride MH. 
Deze technieken zijn een combinatie van twee soorten metaheuristieken. 

We onderscheiden drie types van hybride MH (Blum en Roll, 2003): 

• MH die componenten van een andere MH opnemen 

• cooperatieve zoeksystemen, algoritmen die info uitwisselen 

• integratie van methodes 

Zo zijn er bijvoorbeeld de memetische algoritmen (Corne et al., 1999). Deze zijn 
een combinatie van genetische en local search algoritmen. Meestal manifesteert deze 
combinatie zich als een evolutionair algoritme dat een local search strategie gebruikt 
in plaats van een mutatie-move (Prins, 2004). Bij een evolutionair algoritme wordt 
de oplossingsruimte vrij willekeurig doorzocht. De 'local search-mutatie' zorgt voor 
een meer gerichte zoektocht door de gigantische oplossingsruimte (Rodrigues en 
Ferreira, 2001). 
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Nummer 


Voorwerp 


Gewicht (kg) 


Winstbijdrage 


1 


Tent 


3 


9 


2 


Kleren 


2 


7 


3 


Sandalen 


1.2 


2 


4 


Boek 


0.5 


1 


5 


Kookgerei 


1.1 


9 


6 


Geld 


0.1 


10 


7 


Wandelstokken 


0.6 


5 


8 


Stafkaarten 


0.4 


8 


9 


Slaapzak 


0.75 


5 


10 


Overlevingsgids 


0.9 


10 


11 


matje 


1.2 


4 


12 


snorkel 


0.3 


1 


13 


drinkbus 


1.5 


4 


14 


hamer 


2.5 


2 


15 


stoel 


2.1 


1 



Tabel 2.1: De voorwerpen die de bergbeklininier kan nieenenien, sanien met hun 
gewicht en winstbijdrage 

2.3.3 Voorstelling van een oplossing 



Oni een probleem te kunnen optinialiseren is het noodzakelijk om het op een ge- 
structureerde nianier voor te stellen. Het is interessant om een deel van de struc- 
tuur van een oplossing te behouden. Er zijn enkele fundamentele manieren om een 
optimalisatieprobleem voor te stellen, bijvoorbeeld binaire voorstellingswijze, per- 
mutatievoorstelling. . . 

Bij het eerder besproken klassieke knapsack-probleem is een binaire voorstellingswij- 
ze aanbevolen. Dit kan het best verduidelijkt worden met behulp van een voorbeeld. 

Een bergbeklimmer wil een veertiendaagse trektocht maken door de Al- 
pen. Zijn rugzak heeft een capaciteit van 12kg. Weike voorwerpen moet 
kan hij best meenemen, gegeven het gewicht en de winst van elk van de 
voorwerpen (tabel 2.1)? 

Dit probleem kunnen we voorstellen door een reeks nullen een enen (van- 
daar de naam binair). Een 1 op de ?^-de plaats in de rij betekent dat element 
n wordt meegenomen. Bijvoorbeeld: 
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• In het geval dat er geen enkel object wordt meegenomen: 

000000000000000 

• In het geval dat enkel object nummer 4 en 10 worden meegenomen: 

000100000100000 

Deze voorstellingswijze is duidelijk zeer eenvoudig en toch gestructureerd. Helaas 
zijn niet alle problemen even geschikt oni binair te worden voorgesteld. In het geval 
van het Travehing Salesman Problem is een permutatie-voorstelling aan te raden. 
Hieronder wordt dit verduidelijkt met een voorbeeld. 

Een verkoper van vrachtwagenonderdelen moet dagelijks zes garages (A, 
B, C, D, E en F) bezoeken om ze in te lichten over de nieuwste ontwikkelingen 
op de markt. Gegeven de posities van de steden, wat is de kortste route voor 
de verkoper? 

Een opiossing kan worden voorgesteld door een rij geordende letters of 
getallen. De volgorde van deze letters bepaalt de volgorde waarin de verkoper 
de garages bezoekt. Bijvoorbeeld: 

B A F E C D 

Hierbij begint de verkoper in garage B, om vervolgens langs A, F, E, C 
te gaan en te eindigen bij garage D. 

Er zijn nog een heleboel andere voorstellingswijzen, vb. een matrixvoorstelling die 
binaire voorstellingen combineert. Vaak wordt een representatie van de opiossing 
individueel aangepast aan het probleem. In het uitgewerkt voorbeeld dat we gaan 
bekijken in hoofdstuk 4 zal dit ook het geval zijn. 



2.4 Enkele bekende Metaheuristieken 

2.4.1 Iterated Local Search 

Een metaheuristiek is bij voorkeur zeer simpel, zowel conceptueel als praktisch, 
terwijl hij toch efficient blijft. Iterated local search (ILS) voldoet aan deze eigen- 
schappen (Lourenco et al., 2001). De essentie van ILS is dat het op een iteratieve 
manier oplossingen genereert door een onderliggende heuristische methode. Dit leidt 
tot betere oplossingen dan herhaalde willekeurige runs van die heuristiek. Volgens 
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(MetalieuristicsNetwork, 2004) zijn er twee kernpunten die van een algoritme een 
ILS maken: 

• er wordt een enkele keten gevolgd (dit geldt zeker niet voor populatiegebaseer- 
de algoritnien) 

• de zoektoclit naar betere oplossingen gebeurt in een verkleinde zoekruimte die 
gedefinieerd wordt door de output van een black-box lieuristiek. 

Ondanks zijn zeer simpele aard lieeft liet algoritme een lange gescliiedenis acliter de 
rug. Ondernieer daardoor lieeft liet een reseni verschillende benamingen gekregen, 
zoals iterated descent, large-step Markov-chains, iterated Lin-Kerniglian en chained 
local optimization (Lourenco et al., 2001). 

De working is kort weergegeven in figuur 2.7. We vertrekken van een willekeurige op- 
lossing. Eerst passen we daarop een local search strategie toe zoals die is beschreven 
in sectie 2.3.2. Dit doen we tot we een lokaal optimum bereiken. Dan volgt er een 
'perturbatie', een sprong naar een andere plaats in de zoekruimte. Van daaruit gaat 
de local search strategie verder. Dit proces wordt herhaald tot een vooraf bepaald 
stop-criterium is bereikt (Boston et al., 2001). 

De toegepaste perturbatie moot groot genoeg zijn, dit om to vermijden dat we 
meteen terugkeren naar hot vorige locale optimum. Langs de andere kant mag de 
sprong ook niet te groot zijn (Lourenco et al., 2001). 

2.4.2 Simulated Annealing 

1983 was een belangrijk jaar op hot gebied van combinatorische optimalisatie. Kirch- 
patrick, Gelatt en Vecchi stelden een nieuwe heuristische methode voor, nl. Simu- 
lated Annealing (Gendreau, 2002). Opmerkelijk aan doze techniek is dat men heeft 
aangetoond dat zij steeds converteert naar een optimale oplossing na een oneindig 
aantal moves. Doze ontdekking wekte de interesse van de onderzoekswereld. In de 
jaren die hierop volgden zijn er nog een heleboel verschillende methodes ontwikkeld, 
doze waren ook meestal gebaseerd op natuurkundige processen (Gendreau, 2002). 

Simulated Annealing is gebaseerd op een fysisch proces, nl. hot afkoelen van een 
kristallijne stof. Bij een kristallijne stof zijn alio atomen gerangschikt in een kris- 
talrooster (zie figuur 2.8). Wanneer men de temperatuur doet stijgen, krijgen de 
atomen meer energie en gaan ze sneller bewegen. Hierdoor zitten ze niet meer vast 
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Initiele, willekeurige 
opiossing s 



1 



Willekeurige 
perturbatie s --> t 



I 



Local Search t --> u 
tot u = lokaal optimum 



I 



Als u aanvaard wordt: 
s <-- u 




'Stop- 
criterium 
Neen ^vj^ereikt?- 



Figuur 2.7: Iterated Local Descent algoritme, gebaseerd op (Lourenco et al., 2001) 
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in het rooster. Deze toestand van wanorde kan men vergelijken met de begintoe- 
stand van een optimalisatieprobleem. De structuur van het kristalrooster stelt de 
oplossing van het optimahsatieprobleem voor (Dreo et aL, 2003). 



v\5. 






Figuur 2.8: Na^ en CI m het kristakooster van keukenzout (RUN, 2004) 

Bij het afkoelen van zo'n kristahijne stof kunnen er imperfectheden optreden. Dit 
zijn lokale optima, want elke wijziging zorgt voor een energie-verhoging. De bedoe- 
hng is echter om een globaal optimum te bereiken waarbij de energietoestand zo laag 
mogelijk is. De objectieffunctie komt in dit geval overeen met de mate waarin alle 
atomen op hun juiste plaats zitten in het kristalrooster (^ lage energietoestand). 
Een 'move' is een herschikking van de atomen (daling of stijging van de temperatuur) 
(Kirkpatrick et aL, 1983). 

Het algoritme is voorgesteld in figuur 2.9. We vertrekken van een willekeurige, 
initiele oplossing x, met een begintemperatuur T. We maken vervolgens een random 
move vanuit deze oplossing. De move wordt met zekerheid geaccepteerd als hij de 
oplossing verbetert (lagere energie). Wanneer de oplossing slechter wordt passen we 
de beslissingsregel van Metropolis toe (Kirkpatrick et aL, 1983): de oplossing wordt 
aanvaard met een kans p = e t , 

De kans p om de oplossing te aanvaarden wordt beinvloed door 2 variabelen, de 
temperatuur T en de stijging van de objectieffunctie AE' (Dreo et aL, 2003): 

• Hogere T -^ Hogere kans p (en omgekeerd: [T ^ i p)'- 

AE 

t p = e ^T -L I 

• Hogere Ai? -^ Lagere kans p (en omgekeerd: [T ^ ] p)'- 

I _AET^ 

ip = e T \ I 
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Initiele, willekeurige 
opiossing 



I 



Initiele temperatuur 



1 



Elementaire move 
-> wijziging energie dE 



I 



Aanvaardingsregel van Metropolis: 

- dE = : opiossing aanvaard 

- dE > : opiossing aanvaard met 

kans p = exp (-dE/T) 




Figuur 2.9: Simulated Annealing (Dreo et al., 2003) 
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Een belangrijke observatie is dat men zelf liet aantal aanvaarde moves kan controle- 
ren door T aan te passen. Dit gebeurt door de invoering van een afkoelings schema. 
Enkele mogelijke afkoelingsschema's zijn weergegeven in figuur 2.10. 



Ti=T(To-T,) 
x(1-tanh(^^-5)) 





Figuur 2.10: Mogelijke afkoelingsschema's (Luke, 2004) 



Om de working van dit algoritme beter te begrijpen is hot interessant om simulated 
annealing toe te passen op liet eerder aangehaalde voorbeeld van liet knapsack- 
probleem. 

We gaan verder met het knapsack probleem dat gegeven is in tabel 2.1. 
Stel dat we op een gegeven ogenblik over volgende gegevens beschikken: 

• r = 20 

• de huidige opiossing is: 

011101011110111 

• het gewicht van deze opiossing is 8.2kg. De opiossing is geldig (8.2 < 
12) 

• de winst van deze opiossing bedraagt 54. 

Een mogelijke move die kan worden toegepast is het omdraaien van een 
bit. Stel bijvoorbeeld dat de volgende move het omdraaien van de eerste bit 
is. Dan krijgen we het volgende: 

• de opiossing is: 

111101011110111 
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• het gewicht van deze opiossing is 11.2kg. De opiossing is geldig 

(11.2 < 12) 

• de winst van deze opiossing bedraagt 63. 

Deze move wordt zeker aanvaard omdat hij een daling van de energie teweeg 
brengt {Awinst = AE = 63 — 54 = 9). De temperatuur wordt dus ook 
verlaagd, T = 19 en we gaan verder met de volgende move. Stel dat 
de volgende move het willekeurig omdraaien is van bit 2, dan krijgen we 
volgende opiossing: 

• de opiossing is: 

101101011110111 

• het gewicht van deze opiossing is 9.2kg. De opiossing is geldig (9.2 < 
12) 

• de winst van deze opiossing bedraagt 56. 

Aangezien deze move een stijging van de energie teweegbrengt, wordt 
de beslissingsregel van Metropolis toegepast. De kans p om deze move te 
aanvaarden is: 

p = e~^r~ = e^~w> = 0.623 

Hiervoor trekken we een willekeurig getal u tussen en 1. Als u groter is 
dan p dan zullen we de move verwerpen, en anderszullen we hem aanvaarden. 

Over het proces van Simulated Annealing is reeds veel tlieoretiscli onderzoek ver- 
riclit. Zo is er ondernieer bewezen dat het algoritme na een oneindig aantal moves 
steeds tot de optimale opiossing komt (Blum en Roll, 2003). We moeten ons echter 
wel rekening geven van het feit dat een dergelijk bewijs niet aantoont dat het een 
good algoritme is. Het is in de praktijk belangrijk om snel een kwalitatief goede 
opiossing te vinden, en niet na een oneindig lange tijd de beste opiossing te vinden. 
Dit brengt ons weer naar het belang van metalieuristieken, zoals dit is besproken in 
sectie 2.3. 

2.4.3 Genetische Algoritmen 

De geschiedenis van genetische algoritmen (GAs) begint in 1975, toen Holland het 
book "Adaptation in Natural and Artificial Systems" (Holland, 1975) schreef. Door- 
dat dit book zo'n enorme wetenschappelijke interesse heeft teweeggebracht, wordt 
er wel eens naar verwezen als "The Bible of the GA Community" (Pirlot, 1996). 
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Figuur 2.11: John Holland (Vidal, 2003) 

In tegenstelling tot de hiervoor besproken local search metaheuristieken, zijn GAs 
populatiegebaseerde MH (Pirlot, 1996). Ze behoren tot de klasse van 'evolutionaire 
algoritmen'. Deze evolutionaire algoritnien bevatten naast genetische algoritmen 
ook 'genetische programmering' en 'evolutionary strategies' (Digalakis en Margaritis, 
2001). 

De techniek is sterk geinspireerd op het biologische voortplantingsproces en neemt 
dan ook heel wat ternien uit de biologie over. Hieronder verklaren we de belangrijk- 
ste: 

Populatie Een verzanieling van individuen. 

Individuen Individuele oplossingen, ook wel een chromosomen genoemd (Pirlot, 
1996). 

Fitness Geeft de waarde van de doelfunctie weer voor een individu. Oplossin- 
gen met een hogere fitness zijn beter aangepast aan de onigeving, zij hebben 
nieer kans oni te overleven en hun genen door te geven: "Survival of the fit- 
test" (Blum en Roh, 2003) 

Selectie Het selecteren van individuen uit een populatie die in aanmerking komen 
voor reproductie. Individuen met een hogere fitness zuUen meer kans hebben 
om geselecteerd te worden (Dreo et al., 2003). 

Crossover Het combineren van individuen (parents), om zo hopelijk tot een betere 
individuen (children) te komen. Dit is het reproductieproces (Dreo et al., 2003). 
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Mutatie Een kleine verandering in het genetisch materiaal. Dit is in veel gevallen 
negatief, niaar soms kan het tot een beter individu leiden. Naar dit proces 
refereert men soms als zelfadaptatie (Blum en Roli, 2003). 



Initiele, willekeurige populatie 



Evaluatie van alle individuen 

uit de populatie 

--> fitness 



Evaluatie fffnessvan 
de kinderen 



I 



Selectie van de slechtste 
-> voor vervanging 



Stop 




Figuur 2.12: Genetisch algoritme (Dreo et al., 2003; Towsey et al., 2001) 

Het algoritme is voorgesteld in figuur 2.12. We vertrekken van een initiele populatie 
van n oplossingen. Vervolgens berekenen we voor elk van deze individuen de fitness 
waarde. Op basis hiervan kunnen we de beste oplossingen selecteren met behulp 
van binaire toernooi selectie of roulette wheel selectie. 

Bij roulette wheel selectie, selecteren we een oplossing met een kans die proportioneel 
is aan de fitness- waarde (Dreo et al., 2003). Stel dat we drie individuen in de 
populatie hebben, met respectievelijke fitness waarden: 3, 5 en 6. We gaan deze 
drie oplossingen selecteren met als kans: ^, ^ en ^i- 



In het geval van binaire toernooi selectie selecteren we willekeurig twee individuen 
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en houden we enkel de beste bij. Een belangrijk voordeel hiervan is dat er geen 
risico is om de allerbeste uit de populatie te gooien (Prins, 2004). 

Net zoals het in de natuur ook het geval is, zal bij GAs het meest aangepaste indi- 
vidu overleven en zijn 'genen' doorgeven, daaroni gaan we de geselecteerde (goeie) 
oplossingen combineren. Een volgende stap in het algoritnie is het kruisen van het 
genetisch niateriaal van de geselecteerde individuen (= ouders), via het crossover 
proces. Er zijn verschihende vormen van crossover mogehjk, bijvoorbeeld 1-point 
crossover, 2-point crossover en OX crossover. 

Bij 1-point crossover kiest men 1 kruisingspunt en verwisselen we daarrond het 
genetisch materiaaL Bijvoorbeeld: 



1 1 




1 1 
1 1 1 



10 1111 
10 1 

2-point crossover maakt gebruik van hetzelfde principe, maar dan wel met 2 krui- 
singspunten. Bijvoorbeeld: 

10 110 1 
111 

10 10 1 
1111 

Order crossover (OX) wordt toegepast bij permutatievoorstellingen. Het werkt als 
volgt: 

• We trekken twee kruisingslijnen, net zoals bij 2-point crossover. 

• We behouden het middelste stuk en schrijven dit dus gewoon over van de 
ouders. 

• Voor elk kind: 

— begin na de tweede streep met het invullen van wat er bij de andere ouder 
staat na de tweede streep. 
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controleer voor het invuUen van een cijfer of het niet reeds is gebruikt. 
Dubbele getallen mogen niet voorkomen, in dit geval slagen we het getal 
over en gaan we naar het volgende. 

als we achteraan zijn gekomen, gaan we vooraan verder. 



Bijvoorbeeld: 



1 2 
4 3 



3 4 
2 1 



5 6 

6 5 



2 13 4 6 5 

3 4 2 15 6 

Het resultaat van deze crossover operatie zijn de kinderen. In een volgende stap 
kunnen er willekeurige mutaties optreden. Dit kan bijvoorbeeld het omdraaien van 
een bit zijn. Daarna evalueren we de fitness van de kinderen en selecteren we de 
slechtste individuen uit de populatie, dit kan weerom gebeuren met de niethodes die 
hierboven beschreven staan. Dit proces wordt herhaald totdat een bepaald stop- 
criteriuni voldaan is. 

2.4.4 Tabu Search 




Figuur 2.13: Fred Glover (Glover, 2004) 



Tabu Search (TS) is voor het eerst gepresenteerd door Fred Glover eind jaren '70 
(Dreo et al., 2003). Zijn twee originele artikels (Glover, 1990) (Glover, 1989) bevatten 
veel opvattingen die nog steeds gelden voor tabu search (Dreo et al., 2003). Het 
heeft lang geduurd vooraleer al zijn opmerkingen gebruikt werden in de wetenschap. 
Vooral de aspiratieconditie en het principe van een tabulijst (zie verder in deze 
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sectie) werden al vanaf de jaren '90 gebruikt door wetenschappers die tabu search 
toepasten. Het gebruik van andere structuren liet langer op zich wachten (Glover, 
1998). Op dit moment wordt er nog steeds veel onderzoek verricht naar TS en blijft 
de techniek evolueren en verbeteren (Perry Gray et al., 2002). 

Zoals eerder al bleek is Tabu Search een local search algoritme, met een bijzondere 
eigenschap, het bevat namelijk geheugenstructuren. Glover beschrijft het concept 
tabu search als "a meta-heuristic superimposed on another heuristic" (Perry Gray 
et al., 2002). Het is dus een metaheuristiek die een local search heuristiek helpt 
om een zoekruimte efficient te verkennen. Het basisidee van deze metaheuristiek is 
vermijden dat het zoekproces vast komt te zitten in cirkels. Dit wordt ondermeer 
bereikt door bepaalde moves te verbieden, meer specifiek, diegene die lets doen wat 
juist gedaan is. Deze 'verboden' oplossingen zijn taboe, vandaar de naam tahu search 
(Laguna en Glover, 2004). 

We kunnen een aantal verschillende geheugenstructuren onderscheiden. Fundamen- 
teel kunnen we deze onderverdelen in korte- en lange termijn-geheugenstructuren. 

Korte termijn-geheugenstructuren 

Wellicht de belangrijkste geheugenstructuur is de tabulijst. Deze lijst bevat de laatst 
bezochte oplossingen (dit kunnen eventueel ook moves of attributen van oplossingen 
zijn) en heeft als doel te verhinderen dat we op onze stappen terugkeren (Silver, 
2002). Een eigenschap van de tabulijst is de 'tabu tenure', de lengte van de lijst. 
Wanneer de tabulijst de /c-laatste oplossingen bevat, dan zal de tenure k zijn. De 
moves die we niet mogen uitvoeren omdat de tabulijst ze verbiedt, noemen we 'tabu 
actief. 

Een tweede korte termijn geheugenstructuur is het aspiratieniveau. Hierbij houden 
we steeds de allerbeste oplossing bij. Wanneer we op een gegeven ogenblik tijdens 
het zoekproces een oplossing tegenkomen die beter is dan het aspiratieniveau, dan 
gaan we die zeker aanvaarden (Blum en Roll, 2003). 

Ten slotte bestaan er ook nog de kandidaat lijsten. Deze lijsten bevatten bepaalde 
moves waarvan we op voorhand denken dat ze interessant zijn. De moves die een 
grote bijdrage leveren gaan we vervolgens meer laten voorkomen (Dreo et al., 2003). 
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Lange termijn-geheugenstructuren 

De LT-gelieugens zorgen ervoor dat TS twee belangrijke componenten bezit, nl. 
intensificatie en diversificatie. Tijdens de intensificatie focust de strategie zicli op het 
verkennen van de neighbourhood en de ehte oplossingen. Terwijl bij de diversificatie 
de nadruk hgt op het exploreren van niet verkende ruimten, wat leidt tot oplossingen 
die significant verschillen van diegene die ervoor gekend waren (Blum en Roll, 2003). 

Er zijn vier fundanientele vornien van lange terniijn geheugen. Een eerste type is 
het frequentiegeheugen. Hierbij gaan we bijhouden hoeveel keer een bepaald soort 
oplossing voorkonit. Dit kan bijvoorbeeld worden geiniplementeerd onder de vorm 
van een matrix die bijhoudt hoeveel keer elk attribuut voorkonit. Het kan nuttig 
zijn voor de zoekstrategie om af en toe te starten van een attribuut dat nog niet 
veel is voorgekomen (Dreo et al., 2003). 

Een andere vorm van geheugen is het elitegeheugen. Dit is een lijst met heel goeie op- 
lossingen die we al zijn tegengekomen. Ook hier kan het nuttig zijn om te vertrekken 
van een eliteoplossing als de zoekfunctie vast zit (Sorensen, 2002). 

Vervolgens is er ook de strategische oscillatie. Dit wil zeggen dat we toelaten dat 
de oplossing tijdelijk niet geldig (ook wel 'infeasible' genoemd) wordt. Tenslotte 
hebben we ook de reactive tabu search. Dit is tabu search waarbij de tabu tenure 
kan worden aangepast in de loop van het zoekproces (Wikipedia, 2004). 

De basiswerking van het algoritme is voorgesteld in figuur 2.14. We vertrekken van 
een initiele oplossing s en een lege tabulijst. Vervolgens gaan we de neighbourhood 
genereren. Bij een binaire oplossing kan dit bijvoorbeeld gebeuren door elke bit 
om te draaien. De moves die op de tabulijst staan gaan we editor niet uitvoeren. 
Hierop volgend selecteren we de beste oplossing uit de neighbourhood. Als het 
stop-criterium niet voldaan is gaan we de gedane move op de tabulijst zetten en ver- 
trekkend van de gekozen oplossing de neighbourhood opnieuw genereren, enzoverder 
(Pirlot, 1996). 

Om enkele van deze abstracte begrippen lets duidelijker te niaken, illustreren we het 
tabu search algoritme aan de hand van het in knapsack probleem dat in tabel 2.1 is 
voorgesteld. 

We gaan verder met een knapsack probleem dat gegeven is in tabel 2.1. 
Stel dat we op een gegeven ogenblik over volgende gegevens beschikken: 

• Tabulijst = [5, 7, 8, 9], tenure = 4. 
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Figuur 2.14: Tabu Search Algoritme (Dreo et al., 2003) 
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• de huidige opiossing is: 

011101011110111 

• het gewicht van deze opiossing is 8.2kg. De opiossing is geldig (8.2 < 
12) 

• de winst van deze opiossing bedraagt 54. 

Stel dat het algoritme als move alle bits een voor een gaat omdraaien, 
om zo de neighbourhood te bepalen. De neighbourhood is dan van de orde- 
grootte 0{n). We gaan dan voor alle 15 opiossingen uit de neighbourhood, 
het gewicht en hun winst bepalen. 

De neighbourhood is: 
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We kiezen vervolgens de beste opiossing uit deze 15 (diegene met de 
hoogste winst + feasible) en controleren of deze niet op de tabulijst staat 
(zie rood). Stel bijvoorbeeld dat de beste opiossing werd verkregen door 
het omdraaien van bit 5, dan zouden we deze niet aanvaarden, maar de op 
een na beste opiossing kiezen. Dit is een zeer handig mechanisme om te 
vermijden dat we in cirkeltjes blijven vastzitten. Stel dat de beste opiossing 
die geen tabu is, opiossing 4 is. We nemen deze opiossing over als nieuwe 
opiossing, passen de tabulijst aan en we beginnen terug opnieuw. 
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2.5 Keuze van de juiste metaheuristiek 

Omdat metaheuristieken richtlijnen zijn om een heuristiek te maken, moeten ze 
worden aangepast aan elk probleeni. Bepaalde metaheuristieken zullen ecliter ge- 
niakkelijker aan te passen zijn aan verschillende problenien en sneller tot een goeie 
oplossing konien. Zo zijn er voor het vehicle routing problem een aantal zeer effi- 
ciente tabu search-methodes geschreven en geimplementeerd, maar bijna nog geen 
genetische algoritmen. Prins is de enige die er met een memetisch algoritme (com- 
binatie van genetische algoritmen en local search) de tabu search programma's naar 
de troon steekt (Prins, 2004). 

De keuze van de juiste metaheuristiek is niet altijd even simpel. Er zijn geen goeie 
bestaande theorieen die ons hierbij kunnen helpen (Dreo et al., 2003). Daarom is 
de keuze van een goede MH en de aanpassing van de bijhorende parameters vaak 
afhankelijk van de bekwaamheid en de ervaring van de gebruiker in plaats van de 
exacte toepassing van de voorgeschreven regels. 

Het No Free Lunch-theorema zegt ons bovendien dat er geen methode is die voor 
alle problemen de beste oplossing geeft. 

. . . all algorithms that search for an extremum of a cost function perform 
exactly the same, when averaged over all possible cost functions. In par- 
ticular, if algorithm A outperforms algorithm B on some cost functions, 
then loosely speaking there must exist exactly as many other functions 
where B outperforms A (NoFreeLunchTheorems, 2004). 

Dit wil zeggen dat het niet mogelijk is dat er op een dag een metaheuristiek wordt 
ontwikkeld die de beste oplossingsmethode is voor alle problemen. 
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CAC: enkele projecten 



3.1 Inleiding 



In de vorige twee hoofdstukken hebben we gezien dat we het componeren van mu- 
ziek als een optimaliseringsprobleeni kunnen beschouwen. We hebben ook lets nieer 
gezegd over een bepaald soort oplossingsniethodes voor optinialiseringsproblenien, 
namelijk nietaheuristieken. In dit hoofdstuk gaan we de twee vorige koppelen en 
enkele concrete voorbeelden geven van Computer Assisted Composing-projecten die 
reeds zijn gerealiseerd. De besproken projecten zijn GenJam en SIEDP. We begin- 
nen met een overzicht van enkele meer praktisclie dingen zoals het voorstellen van 
muziek. 



3.2 Hoe gaat CAC in zijn werk? 

3.2.1 Voorstelling van Muziek 

Er zijn verschillende mogelijkheden om muziek voor te stellen met een computer. 
Veruit de bekendste is de 'Common Hierarchical Abstract Representation for Music', 
ofwel de CHA RM -methode. Deze abstracte methode wordt voorgesteld in (Small 
et al., 1990) en is de eerste stap naar een algemeen notationeel systeem voor muziek. 
Ze gaat ervan uit dat elke noot een event is, en dat deze drie eigenschappen heeft: 

Pitch De toonhoogte. 

Time Het tijdstip waarop de noot begin. 

Duration Het aantal tijden dat de noot blijft voortduren. 
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Deze drie eigenscliappen vormen de basiseigenschappen. Ze kunnen nog worden 
uitgebreid met o.a. timbre, intensiteit, amplitude,. . . (Wiggins et al., 1989). Een 
mogelijke voorstelling van een noot met beliulp van liet CHARM -systeem., wordt 
gegeven in figuur 3.1. 
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Figuur 3.1: De Ci/A/^M-voorstellingswijze (Wiggins et al., 1993b) 

De CH A RM --methods is niet gebonden aan een specifieke programmeertaal, en kan 
dus in principe gebruikt worden op elk platform (Reynolds, 1995), van ML (een 
functionele taal) tot Prolog (een logische programmeertaal). 

Een belangrijk voordeel van deze voorstellingswijze is dat ze op een relatief eenvou- 
dige wijze een grammatica-aanpak toelaat (Wiggins en Small, 2000). Grammatica 
doet meteen denken aan een taal en we kunnen ons hierbij dan ook afvragen of 
muziek kan beschouwd worden als een taal. Daarvoor zou er een verband moe- 
ten bestaan tussen liet vlak van de expressie (lioorbare muziek) en liet vlak van 
de inlioud (semantisclie structuur). Dit verband kunnen we ook omsclirijven als 
'betekent'. Het onderzoek omtrent de muzikale code is niet volledig, aangezien er 
niet voldoende begrip is van de vlakken van expressie en inlioud. We kunnen ecliter 
met zekerheid stellen dat er een verschil is tussen gesproken taal en muziek. Langs 
de andere kant is er een sterke neiging om verbanden tussen de twee te zien, zoals 
bijvoorbeeld grammatica (Baroni et al., 1983). 

Een muziekstuk is geen discrete, onafhankelijke verzameling van geordende noten. 
Het bestaat uit 'zinnen', meestal vier maten lang, die op liun beurt liierarcliisch 
geordend zijn. Een muzikale zin wordt gedefinieerd als: "een sequentie die een punt 
bereikt dat een gevoel van pauze of aankomst aangeeft" (Ponsford et al., 1997). We 
moeten met deze zinnen en structuren rekening liouden wanneer we muziek willen 
voorstellen. De C//y4i?M-metliode laat dit toe door een constituent te definieren die 
bestaat uit verscliillende events (Wiggins et al., 1989). Zo'n constituent is dus een 
verzameling noten die samen een zin vormen. 
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3.2.2 De kwaliteit van muziek meten 

Een van de moeilijklieden bij CAC-toepassingen is liet bepalen van de doelfunctie. 
Wat is de kwaliteit van de muziek? Naar dit probleem wordt vaak verwezen als de 
'fitness bottleneck' (Gartland- Jones, 2003). Fitness verwijst naar de benaming van 
de doelfunctie bij genetische algoritnien. Bottleneck duidt op de moeilijkheid om een 
doelfunctie op te berekenen en op te stellen. De bedoeling is dat we een soort van 
score kunnen berekenen of toekennen aan de hand waarvan we de niuzikale kwaliteit 
van liet fragment kunnen kwantificeren (Biles, 1994). 

Er zijn twee fundamenteel verschillende manieren om de 'score' van een muziekstuk 
te bepalen: 

• interactief 

• via een doelfunctie 

De eerste manier is interactief. Individuen gaan de muziek beluisteren en er een 
score aan toekennen. Een belangrijk nadeel liiervan is dat deze score subjectief is 
(Gartland- Jones, 2002). Ze zal bijvoorbeeld sterk afhangen van de gemoedstoestand 
en aandacht (Wiggings et al., 1998). Bovendien gaan de iteraties enorm traag omdat 
liet stuk telkens helemaal beluisterd en beoordeeld moot worden. 

Ten tweede kunnen we aan de hand van formele regels een doelfunctie berekenen. 
Hierbij zijn er geen mensen nodig om de totaalscore te berekenen. Op basis van 
enkele criteria berekent de computer zelf de score van een oplossing. Van deze 
criteria wordt meestal het gewogen gemiddelde berekend op basis van hun belang. 
De totaalscore van deze gemiddelden is onze te optimaliseren waarde (Schoenberger, 
2004). Er dient wel te worden opgemerkt dat de kwaliteit van de vooropgezette 
doelfunctie de resultaten sterk zal beinvloeden. 

In de volgende secties volgen enkele voorbeelden van interactieve en geautomatiseer- 
de doelfuncties. 

We hebben in sectie 1.1 gezien dat een optimaal muziekstuk aan enkele criteria moot 
voldoen zodat het een goede harmonische, melodische en ritmische samenhang heeft 
(Kleeven en Philippi, 2002; Schell, 2002). De moeilijkheid is het bepalen van de 
juiste criteria om deze samenhang te meten. Er zijn vele studies gebeurd, maar er 
bestaat geen algemene doelfunctie. 

De meeste studies zijn gebaseerd op meerstemmige muziek en maken gebruik 
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van criteria op basis van de juiste opeenvolging van akkoorden. De volgende 
maatstaven worden onder andere gebruikt bij liet optimaliseren van meerstem- 
niige niuziek: de juiste akkoordenprogressie, het aantal gebruikte gronddrieklan- 
ken\. . . (Schoenberger, 2004; Bilotta, 2000; B.Pardo, 2002). 

3.2.3 Programmeertalen 

Er zijn een aantal progranieertalen bescliikbaar die speciaal ontworpen zijn oni er 
CAC-programma's mee te schrijven. Enkele voorbeelden: 

• Nyquist 

• C Sound 

• OpenMusic 

We zullen hieronder kort iets vertellen over deze drie programmeertalen. 

Nyquist 

Nyquist is een programmeertaal die speciaal ontwikkeld is voor algoritmisch com- 
poneren en geluidssyntliese. Ze ondersteunt zowel high-level componeren als low- 
level signaalverwerking. Dit is vooral interessant voor componisten die op het 
microscopisch abstractieniveau willen werken (Dannenberg, 2005). Nyquist is 
geimplementeerd bovenop XLisp (een uitbreiding van Lisp). We kunnen Nyquist 
zien als een bibliotheek met Lispfuncties die kunnen worden opgeroepen door func- 
ties in een programma. 

Het grote voordeel van Nyquist is dat componisten de voUedige functionaliteit van 
Lisp kunnen gebruiken. Deze taal wordt veel gebruikt in het onderzoek naar arti- 
ficiele intelli gentle. Anderzijds moeten componisten reeds een goede kennis hebben 
van Lisp vooraleer ze voordeel kunnen halen uit de functionaliteit (Miranda, 2001). 

CSound 

In sectie 1.4.3 is reeds kort iets gezegd over de geschiedenis van CSound. De taal 
CSound is volledig geschreven in C. Wanneer we werken in CSound hebben we 



^ Gronddrieklanken zijn de basisnoten van een toonladder (Geeurickx, 1985). 
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steeds twee bestanden nodig: een 'orchestra file' (.ore) en een 'score file' (.sco). Het 
eerste bestaat uit een verzanieling van instrumenten die de computer vertellen hoe 
hij geluid nioet synthetiseren. Het tweede vertelt de computer wanneer hij het moot 
synthetiseren (Anoniem, 2004). In bijlage B wordt een voorbeeld van een .ore en 
.SCO bestand gegeven. 

OpenMusic 

Gerard Assayag en Carlos Agon van het onderzoekscentrum IRCAM in Parijs, ont- 
wikkelden OpenMusic (Truchet et al., 2003). Het is een systeem voor muziekpro- 
grammering dat het midden houdt tussen een programmeertaal en algemene com- 
positiesoftware. OpenMusic is een visuele interface voor de Lisp taal (CommonLisp 
of CLOS). De taal draait dus bovenop Lisp, net zoals Nyquist (Assayag et al., 1999). 

OpenMusic is objectgeorienteerd. Objecten worden gesymboliseerd door iconen die 
men over het scherm kan slepen. Door het slepen en plakken van deze iconen 
kunnen de meeste operaties gerealiseerd worden (IRCAM, 2005b). Een voorbeeld 
van de visuele programmeeromgeving wordt gegeven in figuur 3.2. 

3.3 GenJam 

3.3.1 Inleiding 

GenJam is ontwikkeld door John A. Biles (Biles, 1994). Het programma imiteert 
een beginnende jazz musicus die solo's leert spelen. GenJam is, zoals de naam reeds 
aangeeft, gebaseerd op een interactief genetisch algoritme (zie sectie 2.4.3). De 
solo's die GenJam produceert komen voort uit toonladders die worden gesuggereerd 
op basis van de gekozen akkoord-progressies. De begeleiding van GenJam bestaat 
uit een standaard ritme en een menselijke mentor die real-time feedback geeft. Uit 
de feedback van de mentor wordt een fitness-waarde berekend. Daarna past GenJan 
verschillende genetische operatoren toe om zijn creatie te verbeteren. In wat hierop 
volgt bespreken we de verschillende onderdelen van GenJam lets meer in detail. 
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Figuur 3.2: De interface van OpenMusic (IRC AM, 2005a) 
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3.3.2 De werking van GenJam 

De 'Genetic Jammer', of kortweg GenJam, is ontwikkeld in een Macintosh/Think 
C omgeving, bovenop een CMU MIDI-toolkit. Een visueel overzicht wordt gegeven 
in figuur 3.3. GenJam leest MIDI sequenties in voor piano, bas en drum. Deze 
instrumenten zuUen zijn solo begeleiden. Verder leest GenJam een 'progressie'- 
bestand in (Biles, 2004). Dit bestand bevat informatie omtrent: 

• ritmisclie stijl 

• aantal solo-refreinen 



akkoordenprogressie 



/Seq 



uenties 



Bas 



Piano 



Drum 




Maat- 
populatie 



Zin- 
populatie 



Figuur 3.3: De arcliitectuur van GenJam (Biles, 2004) 

Op basis hiervan wordt een solo gegenereerd in real-time. De mentor, een persoon 
die hiernaar luistert, zal de verschillende delen van de muziek quoteren met (Biles, 
1996): 

• g: goed 

• b: sleclit 



Op basis hiervan genereert GenJam steeds een verbeterde versie het fragment. 
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3.3.3 Chromosoom-voorstelling 

De manier waarop populaties worden voorgesteld is cruciaal voor het succes van 
een genetisch algoritnie. GenJam gebruikt niet een, maar twee populaties: een van 
niaten en een van 'zinnen'^. Een individu uit de niaat-populatie bestaat uit een 
opeenvolging van MIDI-events, terwijl een individu uit de zin-populatie bestaat uit 
enkele maten. Het grote voordeel hiervan is dat de mogelijkheid bestaat dat een 
bepaalde niaat te herhalen in een zin. Dit geeft de luisteraar vaak het gevoel dat er 
een thema inzit (Biles, 1994). 

In figuur 3.4 wordt een voorbeeld gegeven van de gebruikte voorstelling. Individu 
23 is een willekeurig individu van de zin-populatie. Zoals we zien bestaat deze 
populatie uit vier individuen uit de niaat-populatie, namelijk 57, 57, 11 en 38. In 
functie van het gebruikte genetische algoritnie zullen alle getallen in de figuur binair 
worden voorgesteld. Het geheel van alle individuen in de zin-populatie vormt het 
muziekstuk (Biles, 1996). 



23 



Zin-populatie 




I IVIaat-populatie 

rTjTg 7 5 7 8 7 




Figuur 3.4: Voorbeeldpopulaties (Biles, 1996) 



^Muzikale zinnen zijn opeenvolgingen van nieestal vier niaten, in sectie 4.2 wordt liier meer over 
2;ezegd. 
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3.3.4 Fitness Bottleneck 

Zoals we gezien hebben in sectie 1.2.3 is liet zeer moeilijk om een 'fitness'- of doel- 
functie te bepalen voor een niuzielfstuk. Biles heeft ervoor gekozen oni geen algo- 
ritniische functie hiervoor te voorzien. In plaats daarvan bepaalt een nienselijke 
criticus hoe good het muziekfragment is (Biles, 2004). 

Terwijl de mentor luistert naar de GenJani solo, kan hij ofwel een g (good) ofwel een 
b (sleclit) typen voor elke maat. Elke keer hij een g typt, verhoogt de fitness van de 
betreffende maat en zin. Bij een b verlaagt de fitness. Deze fitness-waarden liggen 
tussen -30 en 30. Op die manier kunnen succesvolle maten de minder succesvoUe 
moeilijk overwelmen (Biles, 1994). 

Hoewel dit een goed systeem is om muziek te beoordelen, zorgt het voor een enorme 
vertraging aangezien de mentor steeds naar het fragment moot luisteren. Hierdoor 
komt er ook een limiet op de grootte van de populaties, het aantal generaties enzo- 
verder. 

3.3.5 Genetische operatoren 

Een genetisch algoritme heeft operatoren nodig om de populatie te laten evolueren. 
In deze korte sectie bespreken we hoe de operatoren van GenJam werken. 

De initiele populatie wordt willekeurig gegenereerd, aan de hand van bepaalde kan- 
sen. Bijvoorbeeld: de kans dat een noot blijft duren is ^ (Biles, 1994). 

De selectie en vervanging gebeuren door een aangepaste vorm van toernooiselectie 
(zie sectie 2.4.3). Er worden willekeurig vier individuen gekozen, deze vormen samen 
een familie. De twee beste individuen van een familie vormen de ouders en de twee 
slechtste worden vervangen door de kinderen die voortkomen uit de combinatie van 
de ouders. Deze combinatie gebeurt door een standaard one-point-crossover (zie 
sectie 2.4.3). Een van de kinderen wordt behouden een het andere kind wordt 
gemuteerd. Enkele mogelijke mutaties op zin- en op maatniveau zijn voorgesteld in 
label 3.1 (Biles, 1994). 
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Mutatie-operator Geniuteerde niaat 

Geen 9 7 5 7 15 15 

Omkering 15 15 7 5 7 9 

Rotatie rechts 15 15 9 7 5 7 

Inverteren (15-waarde) 6 8 15 10 8 15 

Transponeren (vb +3) 12 10 8 10 15 15 



Mutatie-operator 



Geen 

Omkering 

Rotatie reclits (vb.3) 



Geniuteerde zin 



57 57 11 38 
38 11 57 57 
57 11 38 57 



Tabel 3.1: Enkele mutatie-operatoren van GenJam(Biles, 2002) 



3.3.6 Resultaten 

GenJam geeft samen met John A. Biles optredens en lijkt vrij goede muziek te 
produceren. Biles is zich bewust van de fitness bottleneck en heeft een poging 
gedaan om deze te elimineren, door de belioefte aan fitness te elimineren. Dit word 
op de volgende manier gedaan (Biles, 2001): 

• de initiele populatie word 'gevoed' met melodische ideeen uit een database van 
bestaande jazz-fragmenten. 

• er word gebruik gemaakt van een intelligent cross-over operator. 

Het resultaat hiervan is ^Autonomous GenJam'' (Biles, 2001). De resultaten van dit 
nieuwe programma werden gekarakteriseerd als interessanter, meer coherent, men- 
selijker en minder repetitief. Hoewel Autonomous GenJam betere muziek genereert 
is het editor strikt genomen geen genetisch algoritme meer (Biles, 2001), aangezien 
het geen gebruik maakt van fitness. 
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3.4 SIEDP 

3.4.1 Inleiding 

Het 'Intelligent System for Piano Learning Aid', ofwel SIEDP, is een ander soort 
CAC-toepassing. SIEDP componeert geen muziek, maar kan worden gebruikt als 
liulpniiddel bij het componeren. Het systeeni berekent namelijk de beste vingerzet- 
ting oni een niuziekstuk op piano te spelen. Op die nianier kan een coniponist snel 
zien of een gesclireven stuk kan gespeeld worden op een piano. 

De primaire doelstelling van SIEDP is ecliter niet om een hulpniiddel te zijn bij 
het componeren, maar om beginnende pianisten te helpen. Beginnende muzikanten 
verliezen vaak snel hun motivatie omdat ze geen rekening houden met de correc- 
te vingerzettingen. Muziekstukken worden zo moeilijker speelbaar en dit is vaak 
frustrerend (Viana en de Morals Junior, 2003). 

3.4.2 Werking van SIEDP 

SIEDP is ontstaan uit de samenwerking tussen professionele muzikanten, elektri- 
sche ingenieurs en informatici. Het bestaat uit twee hoofdcomponenten (Viana en 
de Morals Junior, 2003): 

• applicatiesoftware die de optimale vingerzetting bepaald. 

• een robotisch hand die de melodie speelt, volgens de gekozen vingerzetting. 

Wanneer we de beste vingerzetting van een melodie willen laten berekenen, volstaat 
het om de melodie in MIDI-vorm door te geven aan SIEDP. Het programma gaat 
vervolgens door middel van een genetisch algoritme (zie sectie 2.4.3) de beste oplos- 
sing proberen vinden. De beste oplossing is de vingerzetting waarbij een minimale 
spierinspanning nodig is. Dit is dus diegene met de kleinste vingerafstand (Viana 
en de Morals Junior, 2003). 

SIEDP werkt met populaties die bestaan uit 50 vingerzettingen. Aangezien het hier 
gaat om een genetisch algoritme, is het ook hier weer nodig om een fitness-functie 
te definieren. Deze functie is een maat voor de 'afstand' van de vingerzetting van 
een gegeven melodie. Op basis van deze waarde gaat het algoritme goede individuen 
selecteren uit de populatie. Door selectie, crossover en mutatie gaat de populatie 
evolueren en zo hopelijk tot een goede oplossing komen. 
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3.4.3 Resultaten 

Als output geeft SIEDP een opsomming met de tien beste vingerzettingen. De 
eindgebruiker kan dan zelf nog kiezen welke hij wil gebruiken. In de vorige versie 
van SIEDP was er een probleeni met de convergentie. Het duurde lang voordat 
men een goede oplossing vond. Dit probleem ontstond doordat er veel 'klonen' in 
de populatie zaten, of met andere woorden, bepaalde oplossingen kwamen meer dan 
een keer voor in de populatie. In de laatste nieuwe versie is dit opgelost door een 
routine in te voeren die klonen verwijdert. Hierdoor komt men sneller tot een goede 
oplossing (Viana en de Morals Junior, 2003). 



DEEL II 



Toepassing: Tabu Search voor de optimalisatie 
van muzikale fragmenten 



T" Een tabu search algoritme 
voor het optimaliseren 
van muzikale fragmenten 



4.1 Inleiding 

. . . les calculs n'ont cesse d'accompagner la musique tout au long de son 
histoire, et dans toutes les civilisations . . . 

(Boulez, 1985), zoals aangeliaald in (Harris et al., 1993) 

We hebben zelf een eenvoudige CAC-toepassing geschreven die in staat is muziek- 
fragnienten te optimaliseren. In dit deel bespreken we het zelf ontworpen Pascal- 
programma 'TabuMusic', waarvan de broncode bescliikbaar is in bijlage D en op de 
cd-rom. Op de bijliorende cd-rom staat tevens een geconipileerde versie voor Win32 
en linux. 0ns hoofddoel bij het niaken van TabuMusic was een algoritme te ontwer- 
pen dat een willekeurig gegenereerd muziekstuk optimaliseert met behulp van een 
tabu search-algoritme. Hierbij wordt er een score toegekend aan de muziek. Ver- 
volgens wordt de enorme zoekruimte van mogelijke noten en ritmes op een gerichte 
manier verkend, om zo hopelijk tot een goede oplossing te komen. 

In de bestaande literatuur omtrent Computer Assisted Composing bestaan er weinig 
of geen toepassingen die gebruik maken van tabu search. Het leek ons daarom 
interessant om na te gaan of een tabu search algoritme ook kan worden gebruikt 
om een muziekstuk te componeren. Zoals in hoofdstuk 3 bleek, worden andere 
metaheuristieken, zoals bijvoorbeeld genetische algoritmen, reeds veel gebruikt in 
dit vakgebied. Onze bedoeling is om na te gaan of een tabu search algoritme in 
staat is om snel een good muziekstuk te vinden, gegeven enkele beperkingen. 
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We gaan liomofone muziek optimaliseren, dit wil zeggen dat in het te optimaliseren 
muziekstuk steeds maar een noot tegelijk klinkt. Deze vereenvoudiging zorgt on- 
der meer voor een gemakkelijkere representatie en een sinipelere doelfunctie, zoals 
duidelijk wordt secties 4.2 en 4.4. 

De structuur van de ontwikkelde toepassing is voorgesteld in figuur 4.1. We vertrek- 
ken van een initieel willekeurig muziekfragment en een lege tabulijst. Vervolgens 
kunnen we steeds kiezen of we het ritnie of de melodie willen optimaliseren. De 
gebruiker specificeert op voorhand hoeveel keer hij het ritnie en de melodie wil opti- 
maliseren en in welke volgorde. In wat volgt worden de verschillende stappen meer 
in detail uitgelegd. 



Initieel, willekeurig 
muziekfragment 




Figuur 4.1: Schematische voorstelling van de ontwikkelde toepassing 
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4.2 Voorstellen van muziek 

Een eerste stap is het weergeven van muziek. Hoewel de CHARM -iiotatie (zie sectie 
3.2.1) een goede algeniene voorstellingswijze vormt hebben we er niet voor gekozen 
oni ze te gebruiken. Voor onze toepassing is het belangrijk oni muziek op een simpele 
manier voor te stellen. Een eenvoudige structuur is hiervoor beter geschikt, zodat 
we een duidelijk en gestructureerd beeld krijgen van wat er op elk tijdstip gebeurd 
in het muziekstuk. Aangezien we geen gebruik maken van ingewikkelde regels in 
verband met grammatica, zou de CHARM -m.eth.ode onnodige complexiteit aan het 
model toevoegen. 

We hebben zelf een simpele matrixvoorstelling ontworpen om muziek weer te geven. 
Een muziekstuk wordt hierbij voorgesteld als een eenvoudige (2 x ?72)-matrix. Waar- 
bij 771 de lengte van het fragment is, uitgedrukt in 16de noten. De noot in kolom j 
hoort bij het j-de tijdstip. 

In de eerste rij wordt elke noot voorgesteld door een integer waarde. Deze waarde 
is de midi-waarde van de noot. Een bepaalde toonhoogte wordt voorgesteld door een 
integer, zo krijgt de lage do de waarde 60 toegewezen. Een verhoging van een halve 
toon staat gelijk aan de vorige midi-waarde + 1 (Miranda, 2001). De oorspronkelijke 
bedoeling van deze voorstellingswijze was dat het voor een gemakkelijkere conversie 
naar het MIDI-formaat zou zorgen. Achteraf is deze conversie echter gebeurd via de 
ABC-notatie (zie sectie 4.7.3). 

De tweede rij bevat steeds een 1 of 0, al naar gelang de noot wel of niet blijft voort- 
duren op het volgende tijdstip. Een voorbeeld van de ontwikkelde matrixvoorstelling 
wordt gegeven in figuur 4.4. 

Om realistisch te blijven moot er ook een beperking komen met betrekking tot het 
aantal mogelijke toonhoogtes. Het toonbereik van het muziekstuk gaat van de lage 
do tot en met de hoge si, en omvat twee voUedige octaven. Zo kunnen we ongeveer 
evenveel klanken weergeven als een klassieke fiuit (hoewel deze net nog iets hoger 
kan gaan). De volledige range samen met de overeenkomstige midi-waarden zijn 
voorgesteld in figuur 4.2. 

Een volgende beperking die nodig is, heeft betrekking op het ritme. Tijdens een 
muziekstuk houdt men het ritme aan door te 'tellen'. De duur van de noten is 
gebaseerd op deze tellen en op de maatstaf. Wij hebben als maatstaf | genomen, 
dit wil zeggen dat de totale lengte van een maat bestaat uit 4 noten met duur 
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Figuur 4.2: Gebruikte noten en hun MIDI-waarde 

I (nootduur ^ konit dus overeen met 1 tel). Deze noten kunnen natuurlijk nog 
verder worden opgesplitst in vb. 2 aclitste-noten (| + | = \)- De duur van 1 tel (1 
vierdenoot) is relatief en afhankelijk van het gekozen tempo (dit kan vb. zijn ^ = 120 
keer per minuut). Er bestaan verschillende soorten noten die elk een andere duur 
hebben: 

• Een aclitste-noot | 

• Een kwartnoot ^ 

• Een halve noot ^ 



De basis-nootduren zijn voorgesteld in figuur 4.3. Een basisnootduur kan aangepast 
worden door er een puntje naast te zetten, dan doet men zijn duurxl.5. Er zijn 
nog andere aanpassingen mogelijk om nog meer verschillende nootduren te kunnen 
voorstellen. 
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8 


4 


16 


8 


16 



Figuur 4.3: Ritmedeeltjes 

Om de voorstelling van het te optimaliseren stuk te vereenvoudigen nemen we als 
kleinst mogelijke ritmedeeltje de 16e noot en als langste de hole noot. Samenvattend 
kunnen we stellen dat een muziekstuk van n maten zal bestaan uit maximaal n x 16 
noten. Het voUedige muziekstuk komt in een 2 x (n x 16)-matrix. Een voorbeeld 
hiervan is uitgewerkt in figuur 4.4, waar een kort fragment van Senate KV 545 van 
Mozart in matrixvorm is beschreven. 



In deze matrix hebben de elementen volgende betekenis: 
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i 



P 



A = 



72 72 72 72 72 72 72 72 76 76 76 76 79 79 79 79 
1111111011101110 
71 71 71 71 71 71 72 73 71 71 71 71 -1 -1 -1 -1 
1111100011101 1 1 



Figuur 4.4: Wolfgang Amadeus Mozart - Sonate KV 545 

• Element ai^ van de matrix geeft de midi-waarde van de noot op de k-de positie. 
Wanneer het om een rust gaat is deze waarde —1. 

• Element a2k van de matrix geeft aan of de /c-de noot blijft voortduren (1) of 
niet (0). Een duidt erop dat er ofwel een andere noot/rust volgt ofwel dat de 
noot herhaald wordt (en dus niet blijft voortduren). Elke kolom staat voor een 
7^-noot. Dit is dus de kleinste ritmische eenheid die we kunnen voorstellen. 

16 



4.3 Initiele willekeurige muziek 



We vertrekken van een willekeurig muziekfragment, dat we later gaan optimaliseren. 
Hoe genereren we deze muziek? Dit zal gebeuren op basis van een aantal parameters 
die de gebruiker doorgeeft aan het systeem: 

• de toonaard. 

• het aantal maten muziek dat hij wil genereren. 

Wat betreft de toonaard kunnen we opmerken dat alle mogelijke majeur-toonaarden 
mogelijk zijn. Een uitbreiding naar mineur of andere speciale toonaarden kan later 
gemakkelijk worden gerealiseerd. Daarvoor hoeven we enkel de manier waarop de 
Hoonladderverzameling^ wordt gegenereerd aan te passen. Deze toonladderverzame- 
ling is een simpele matrix die alle mogelijke noten bevat uit de toonladder. Deze 
matrix wordt helemaal in het begin van het programma gegenereerd. 

De parameter ^aantal maten muziek die de gebruiker wil genereren', spreekt voor 
zich. De gebruiker kiest op voorhand de lengte van het muziekfragment. Dit blijft 
constant gedurende het optimalisatieproces. 
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Het genereren van een willekeurig muziekstuk gebeurt zoals is voorgesteld in fi- 
guur 4.5. Voor elke noot wordt er eerst een toonhoogte gekozen. Een rust 
wordt beschouwd als een noot met toonhoogte -1, deze is niet opgenonien in 
de toonladderverzameling. De kans p oni een bepaalde toonhoogte te kiezen is 
p = l/(aantal noten in toonladderverzamehng + 1). Vervolgens krijgt elke noot een 
duur. We hebben de keuze uit 16 mogelijke nootduren (zie sectie 4.2). De kans q op 
een bepaalde nootduur is dus q = j^. 



Voor elke noot: 



I 



1. Random toonhoogte kiezen: 



Kans p = (a) / {a+1) 

--> toonhoogte uit toonladderverzameling 
Kans p = 1 / {a+1) 

--> rust (= noot met toonhoogte -1) 



T 



2. Random nootduur k\ezer\: 



Kans p = 1 / 1 6 

--> een van de 16 mogelijke nootduren 



Neen 




Noot afbreken indien nodig 

-> aanpassen aan lengte 

muziekfragment 



I 



Stop 



a = aantal noten in de toonladderverzameling 



Figuur 4.5: Algoritme generatie willekeurig muziekfragment 



Nadat het willekeurig muziekstuk is gegenereerd zal er een score moeten toegekend 
worden aan het fragment. Deze geeft de muzikale kwaliteit weer van de gegenereerde 
muziek. 
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4.4 Muzikale kwaliteit van een fragment 

We hebben gekozen voor een geautomatiseerd doelfunctie, omdat deze de meeste 
voordelen biedt. Een doelfunctie berekenen is zeer snel en geeft een objectief resul- 
taat. 

Als criteria voor het nieten van de kwaliteit van lioniofone niuziek vinden we bij 
verscliillende auteurs (Snarrenberg, 2004; Scliell, 2002; Papadopoulos en Wiggins, 
1998) maatstaven terug die sterk met elkaar verband liouden o.a. geconnecteerdheid, 
varieteit, intervalgrootte, consonantie, toonbereik,. . . Wij gebruiken de criteria van 
Towsey omdat deze de meeste van de vorige omvatten (Towsey et al., 2001). 

De gebruikte criteria kunnen elk worden voorgesteld in breukvorm: (Towsey et al., 
2001): 

Pitch variety De verhouding tussen het aantal verschillende toonhoogtes en het 
aantal noten. Dit criterium geeft de diversiteit in het aantal klanken weer. 

•J • 1 aantal verschillende toonhoogtes 

criterium 1 = --; — : — 

aantal noten 

Dissonant intervals De fractie dissonante intervallen. De dissonantiewaarde van 
de verschillende intervallen (gemeten in halve tonen en de rusten genegeerd), 
worden gegeven in tabel 4.1. 



Interval 


Dissonantiewaarde 


0, 1, 2, 3, 4, 5, 7, 8, 9, 12 


0.0 


10 


0.5 


6, 11, > 13 


1.0 



Tabel 4.1: Dissonantiewaarden 



criterium 2 = 



^ alle dissonantiewaarden 
aantal intervallen(=aantal noten — 1) 



Contour direction De verhouding tussen het aantal stijgende intervallen en het 
totaal aantal intervallen. Rusten worden genegeerd. 



• j_ • o y^. alle stii 

criterium o = — — v^— it 



rende intervallen 



^ alle intervallen 
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Criterium 


Gemiddelde 


Standaardafwijking 


Pitch variety 


0.27 


0.11 


Dissonant intervals 


0.01 


0.02 


Contour direction 


0.49 


0.06 


Contour stability 


0.40 


0.11 


Rhythmic variety 


0.24 


0.07 


Rhythmic range 


0.32 


0.11 



Tabel 4.2: De belangrijkste criteria samen met hun gemiddelde en standaardafwij- 
king 



Contour stability De proportie intervallen die in dezelfde richting evolueren als 
het vorige interval. Rusten worden genegeerd. 

-J • A aantal opeenvolgende intervallen die in dezelfde richting evolueren 



aantal intervallen — 1 



Rhythmic variety De ritmische coherentie: de graad van gebruik van de 16 mo- 
gelijke ritmes. 

aantal verschillende nootduren 



criterium 5 = 



16 



Rhythmic range Het ritmische bereik 



criterium 6 = 



max nootduur — min nootduur 
16 



Deze criteria alleen zijn echter niet voldoende. We moeten ook hun optimale waarde 
weten. Towsey heeft hiernaar een onderzoek gedaan. Hij heeft een aantal bekende 
muziekstukken uit verschillende muziekstijlen ingevoerd in de computer en de waarde 
berekend voor elk bovenstaand criterium en een aantal andere criteria (Towsey et al., 
2001). De belangrijkste hiervan (die met de kleinste standaardafwijking) zijn diegene 
die net besproken zijn geweest. Ze staan opgesomd in tabel 4.2, samen met hun 
gemiddelde waarde en standaardafwijking. 

We vormen de criteria om zodat ze minimaal zijn wanneer hun optimale waarde (het 
gemiddelde in tabel 4.2) bereikt wordt. Om de doelfunctie concreet te berekenen 
nemen we een gewogen gemiddelde. De gebruiker kan de wegingscoefficenten als 
parameter doorgeven aan het programma. In sectie 5.2 bespreken we de keuze van 
deze wegingscoefficenten. We krijgen dan volgende (gestandaardiseerde) functie: 
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|criteriuni 1 — 0. 27| . |criterium 2— 0.01| . |criterium 3— 0.49| . |criteriuiTi 4 — 0.40| . 

score — tti Q-^j^ h a2 q^qj r fls g^oe r a^ ^^-^ h 

criterium 5— 0.24| j^ |criterium 6— 0.32| 

^5 om I" ^6 0:11 

met tti = wegingscoefficient van criteriuni i 

Deze functie is optimaal wanneer ze nul is. Daarom is het doel van ons tabu search- 
algoritnie oni deze functie te minimaliseren. 



4.5 Local search operators 

We optimaliseren de willekeurig gegenereerde muziek met behulp van een tabu search 
algoritme. Zoals in sectie 2.4.4 is uiteengezet, hebben we steeds een local search- 
operator nodig om onze neighbourhood te kunnen berekenen. Enkele voorbeelden 
van mogelijke moves kwamen reeds aan bod in sectie 3.3. In wat volgt bespreken 
we de local search-operators van de melodie en het ritme. 

4.5.1 Optimalisatie melodie 

Als move kozen we het vervangen van een noot door een noot met dezelfde leng- 
te, maar met een andere toonhoogte. De neighbourhood van een muziekfragment 
bestaat uit alle fragmenten waarbij telkens een noot wordt vervangen. Deze noot 
wordt in de verschillende fragmenten vervangen door alle mogelijke noten uit de 
gekozen toonladderverzameling. We proberen alle noten te vervangen. 

Een voorbeeld maakt dit lets concreter. 

Stel je wil volgend (zeer) kort fragment optimaliseren: 



fe 



Het is een fragment in C, wat wil zeggen dat er geen tj of b voorkomen 
in de toonladderverzameling. De toonladderverzameling is hieronder voor- 
gesteld. We merken op dat ze alle noten bevat die mogen voorkomen in het 
fragment en die tot de toonladder behoren. 

^ o o o " ' = 
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Hieronder wordt de hele neighbourhood voorgesteld. Een move wordt 
weergegeven als een maat in de figuur. 






li. J J 1 - J 1 ,. J 1 F J 1 ^^ 


J^ — ^ si si si si si 



/p iI p iI p iI p II p^ 

g r j ' r J ' r J ' r -J ' ^ 



L ? r. \ ? P I P P I P [^ I ^^ 



4.5.2 Optimalisatie ritme 

Een move is het veranderen van de duur van een noot door eender welke andere 
niogelijke duur. De neighbourhood van een niuziekfragnient bestaat dus uit alle 
fragmenten waarbij telkens de lengte van een noot wordt aangepast. We proberen 
alle noten aan te passen. 

Een voorbeeld maakt dit lets concreter. 

We optimaliseren een zeer kort fragment dat bestaat uit drie noten, zie 
tabel 4.3. Om te verduidelijken hoe we een move maken stellen we het 
fragment in tabelvorm voor: 



Iste noot duur 



sol 



_8_ 
16 



2de noot duur 



la 



2 
16 



3de noot duur 



do 



_6_ 
16 



Tabel 4.3: Voorbeeldfragment 

De nootduur wordt steeds voorgesteld als een veelvoud van j^. Dit komt 
omdat we ^ als kleinste ritmedeeltje hebben beschouwd (zie sectie 4.2). 
De som van alle nootduren in een maat moet steeds y| zijn. Aangezien 
ons voorbeeldfragment maar uit een maat bestaat moet de som van alle 
nootduren steeds y| zijn. Op deze manier blijft de lengte van de muziek 
steeds even lang. 
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Een move is het aanpassen de nootduur van een noot. Voor eike noot 
zullen we alle mogelijke duren proberen. De neighbourhood die hierbij hoort 
is voorgesteld in tabel 4.4. 

We zien echter dat de som van alle nootduren groter is dan y|. Er zijn 
twee mogelijke gevallen: 

1. lengte > jf ^ kap de laatste noot af. 

2. lengte < ^f ^ vul het muziekstuk achteraan aan met de verkorte noot. 
Dit is uitgewerkt in tabel 4.5. 



4.6 Tabu-structuren 

De algoritmen voor het optimaliseren van het ritme en de melodie hjken sterk op 
elkaar. In figuur 4.6 worden ze scheniatisch weergegeven. In beide gevahen gaat het 
om een tabu search algoritnie. 

4.6.1 Optimalisatie melodie 

De optiniahsatie van de melodie gebeurt door het oproepen van de functie noot- 
transponeren. De working van deze tabu search functie is scheniatisch voorgesteld 
in figuur 4.6. 

De allereerste keer vertrekken we van het willekeurig gegenereerde muziekstuk en 
een lege tabulijst. De gebruiker kan in het begin volgende parameters in verband 
met de melodie doorgeven aan het systeem: 

• de lengte van de tabulijst, ofwel de tabu tenure. 

• het aantal iteraties, hoeveel keer willen we het algoritnie doorlopen en een 
move niaken? 

Eerst wordt de neighbourhood van het originele fragment s gegenereerd zoals is 
uiteengezet in de vorige sectie. Vervolgens wordt de muzikale kwaliteit van alle 
oplossingen uit de neighbourhood berekend. Hoe good is hun score? Wanneer een 
oplossing beter of even good is dan de score van s, dan wordt dit onze nieuwe 
oplossing t. Dit zal echter enkel gebeuren als de move niet op de tabulijst staat. In 
de tabulijst houden we bij welke noten onlangs vervangen zijn en welke midi-waarde 
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duur 111 


duur n2 


duur n3 
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Tabel 4.4: Voorbeeld neighbourhood 
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duur 111 


duur n2 


duur n3 


nieuwe n4 


duur 111 


duur n2 


duur n3 


nieuwe n4 


1 

16 


2 
16 


6 
16 


Yg noot 1 


8 
16 


8 
16 


/ 




2 

16 


2 
16 


6 
16 


yI noot 1 


8 
16 


8 
16 


/ 




3 
16 


2 
16 


6 
16 


^ noot 1 


8 
16 


8 
16 


/ 




4 
16 


2 
16 


6 
16 


il noot 1 


8 
16 


8 
16 


/ 




5 
16 


2 
16 


6 
16 


re noot 1 


8 
16 


8 
16 


/ 




6 
16 


2 
16 


6 
16 


I6 i^oot 1 


8 
16 


8 
16 


/ 




7 
16 


2 
16 


6 
16 


I6 i^oot 1 


8 
16 


8 
16 


/ 




8 
16 


2 
16 


6 
16 




8 
16 


8 
16 


/ 




9 

16 


2 

16 


5 
16 




8 
16 


2 
16 


1 

16 


^ noot 3 


10 
16 


2 

16 


4 
16 




8 
16 


2 
16 


2 
16 


^ noot 3 


11 
16 


2 
16 


3 

16 




8 
16 


2 
16 


3 
16 


^ noot 3 


12 
16 


2 
16 


2 
16 




8 
16 


2 
16 


4 
16 


^ noot 3 


13 
16 


2 
16 


1 

16 




8 
16 


2 
16 


5 
16 


j^ noot 3 


14 
16 


2 
16 


/ 




8 
16 


2 
16 


6 

16 




15 
16 


1 
16 


/ 




8 
16 


2 
16 


6 
16 




16 
16 


/ 


/ 




8 
16 


2 
16 


6 
16 




8 
16 


1 
16 


6 
16 


Yq noot 2 


8 
16 


2 
16 


6 
16 




8 
16 


2 
16 


6 
16 




8 
16 


2 
16 


6 
16 




8 
16 


3 

16 


5 
16 




8 
16 


2 
16 


6 
16 




8 
16 


4 
16 


4 
16 




8 
16 


2 
16 


6 
16 




8 
16 


5 
16 


3 

16 




8 
16 


2 
16 


6 
16 




8 
16 


6 
16 


2 
16 




8 
16 


2 
16 


6 
16 




8 
16 


7 
16 


1 

16 




8 
16 


2 
16 


6 
16 




8 
16 


8 
16 


/ 




8 
16 


2 
16 


6 

16 





Tabel 4.5: Voorbeeld aangepaste neighbourhood 
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ze hadden. Stel dat we op een bepaald ogenblik de volgende tabulijst hebben met 
als tenure 5: 



A = 



173 5 
65 60 



In liet voorbeeld werden reeds twee noten vervangen. De eerste noot die vervangen 
werd begon op de vijfde tijd (vijfde positie in de ontworpen niatrixvoorstelling). 
Deze noot had als midi-waarde 60, het was dus een do. De volgende vervangen noot 
had als midi-waarde 65 en begon op de 173ste tijd. 

In TabuMusic is ook een aspiratieconditie ingebouwd. We houden steeds de allerbes- 
te oplossing bij die we gevonden hebben. Wanneer we een oplossing u tegenkomen 
die beter is dan de allerbeste, dan nemen we deze steeds aan. Zelfs als de move op 
de tabulijst staat wordt hij toch gedaan omdat we onze allerbeste oplossing kunnen 
verbeteren. 

Wanneer er een oplossing t oi u wordt aangenomen, dan komt de gedane move op 
de eerste plaats van de tabulijst te staan en schuift alles een plaats op. In het geval 
dat er geen betere oplossing gevonden is, stijgen we naar de minst slechte oplossing. 

Om de volgende move te maken vertrekken we van de laatst aangenomen oplossing. 
Hiervan berekenen we weer de neighbourhood, enzoverder. De gebruiker kan in het 
begin specifieren hoeveel keer de procedure herhaald wordt. 

4.6.2 Optimalisatie ritme 

De optimalisatie van het ritme gebeurt door het oproepen van de functie ritmeverhe- 
teren. De working van deze tabu search functie is schematisch voorgesteld in figuur 
4.6. De gebruiker geeft weer dezelfde parameters door als bij de optimalisatie van 
de melodie. 

Net zoals bij de melodie vertrekken we van het willekeurig gegenereerde muziekstuk 
en een lege tabulijst. Dan wordt de neighbourhood van het originele fragment s 
gegenereerd zoals is uiteengezet in sectie 4.5. Vervolgens wordt de muzikale kwa- 
liteit berekend van alle oplossingen uit de neighbourhood. Hoe good is hun score? 
Wanneer een oplossing beter of even good is dan de score van s, dan wordt dit onze 
nieuwe oplossing t. Dit zal, net zoals bij de melodie, enkel gebeuren als de move niet 
op de tabulijst staat. Wanneer een oplossing beter of even good is dan de vorige. 
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Muziekfragment s 



i 



Bereken de neighbourhood van s 




Figuur 4.6: Optimalisatie van de melodie en het ritme 
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dan zal ze worden aangenonien als de move niet op de tabulijst staat. In de tabu- 
lijst liouden we bij welke noten onlangs vervangen zijn. Stel dat we op een bepaald 
ogenblik de volgende tabulijst hebben met als tenure 5: 



A = 



10 8 2 



Er werden reeds drie noten vervangen. De noot die eerst vervangen werd was de 
tweede noot. Let op, dit is niet de noot die op de tweede tijd begint, zoals bij de 
melodie. Dit is de tweede noot van liet muziekstuk. Vervolgens werd ook de aclitste 
en de tiende noot vervangen. 

Ook voor de optimalisatie van liet ritme is er een aspiratieconditie ingebouwd. We 
liouden steeds de allerbeste oplossing bij die we gevonden hebben. Wanneer we een 
oplossing u tegenkomen die beter is dan de allerbeste, dan nemen we deze steeds 
aan. Zelfs als de move op de tabulijst staat wordt liij toch gedaan omdat we onze 
allerbeste oplossing kunnen verbeteren. 

Wanneer er een oplossing t oi u wordt aangenomen, dan komt de gedane move op 
de eerste plaats van de tabulijst te staan en schuift alles een plaats op. In liet geval 
dat er geen betere oplossing gevonden is, stijgen we naar de minst sleclite oplossing. 

Om de volgende move te maken vertrekken we van de laatst aangenomen oplossing. 
Hiervan berekenen we weer de neighbourhood, enzoverder. De gebruiker kan in het 
begin specifieren hoeveel keer de procedure herhaald wordt. 

4.7 Implementatiedetails 

In deze sectie wordt dieper ingegaan op de specifieke implementatiedetails van Ta- 
buMusic. Welke parameters kunnen we doorgeven, in welke taal is het geschreven 
en hoe kunnen we de gecomponeerde niuziek beluisteren? 

4.7.1 Parameters 

Wanneer we het programma 'TabuMusic' oproepen zal het eerst het bestand 'input- 
muziek.txt' aanroepen. Hierin staan volgende parameters gespecificeerd: 
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Integer 


Toonladder 


60 


C 


61 


ctt 


62 


D 


63 


Dtt 


64 


E 


65 


F 


66 


n 


67 


G 


68 


Gtt 


69 


A 


70 


Att 


71 


B 



Tabel 4.6: De integer-waarden van de verschillende toonladders 



• De toonladder van de te genereren muziek. Dit wordt aangegeven door een 
integer. Tabel 4.6 geeft een overziclit van de mogelijke waarden. Het gaat liier 
steeds om de majeur toonladder. 

• Het aantal niaten muziek dat de gebruiker wil genereren. 

• De bestandsnaam waarin de gebruiker de willekeurige muziek wil opslaan. Dit 
is bij voorkeur een .abc-bestand. 

• De bestandsnaam waarin de gebruiker de geoptimaliseerde muziek wil opslaan. 
Dit is bij voorkeur een .abc-bestand. 

• De bestandsnaam waarin de tekstoutput van TabuMusic komt. 

• Het aantal keer dat de gebruiker wil itereren per optimalisatie. Dit is het 
aantal keer dat er een neighbourhood zal worden gegenereerd en onderzocht 
per tabu search. 

• De lengte van de tabulijst van de melodie, maximaal 1000. 

• De lengte van de tabulijst van het ritme, maximaal 1000. 

• De zes wegingscoefficienten. een voor elk criterium van de doelfunctie. 

• De sequentie van de optimalisaties: dit is een opeenvolging van ofwel de me- 
lodie (1) of het ritme (2) en kan er bijvoorbeeld als volgt uitzien: 12 2 2 11 
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2 1 2. 

• Het bestand wordt afgesloten met een 0. 

TabuMusic leest al deze parameters in wanneer het wordt aangeroepen. Een voor- 
beeld van een inputmuziek-file is bijgevoegd in bijlage C. 

4.7.2 Programmeertaal 

In sectie 3.2.3 zijn enkele talen genoemd die speciaal ontworpen zijn voor CAC- 
toepassingen. Deze talen zijn echter zeer ingewikkeld, daarom hebben wij er geen 
gebruik van gemaakt. Wij kozen ervoor om de gestructureerde taal Pascal te ge- 
bruiken omdat zij volgende voordelen bood: 

• het is een zeer snelle taal 

• ze is zeer gestructureerd 

• de syntax is relatief eenvoudig 

Deze eigenschappen zorgen ervoor dat we de gestructureerde en snelle aard van onze 
oplossingsmethode (tabu search) niet verliezen in een warboel aan programeerob- 
jecten. Bij Pascal zijn er echter geen ingebouwde functies om muziek om te zetten 
in geluid. In de volgende sectie wordt verklaard hoe we dit hebben opgelost. 

4.7.3 De muziek beluisteren 

Wanneer TabuMusic twee muziekstukken gecomponeerd heeft, een willekeurig en 
een geoptimaliseerd, is het interessant om ze aan een menselijk oor te onderwerpen. 
We willen weten of er een verbetering te horen is in het muziekstuk. 

Er zijn verschillende bestandsformaten voor muziek: mp3, wav, midi,...Voor ons 
brengt het midi formaat de meeste voordelen mee. In sectie 1.4.3 word de MIDI- 
standaard reeds meer in detail uitgediept. Voor ons is vooral van belang dat er 
zoiets bestaat als een MIDI-bestand. Dit is afspeelbaar op gelijk welke MIDI-speler. 

Praktisch is niet gemakkelijk om een tekstbestand om te zetten naar een MIDI of 
.mid bestand. Er bestaan programma's die dit kunnen doen voor ons, zie (de Bruijn, 
2004). We hebben echter gekozen om de muziek via een omweg om te zetten in 
midi, door gebruik te maken van de ABC-notatie. Deze notatie stelt elke noot 
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voor door een letter van A tot G (in kleine letters voor het 2e octaaf). Een rust 
wordt voorgesteld als een z. Deze letter wordt steeds gevolgd door een integer die 
de lengte aangeeft. Dit getal is steeds een veelvoud van de basisduur (bij ons ^q) 
(Mansfield, 2004). In figuur 4.7 wordt een niuziekfragnient in matrixvorm omgezet 
in de ABC-notatie. 
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Figuur 4.7: Matrixvoorstelling omgezet in ABC-notatie 

In de header van het ABC-bestand staat nog additionele informatie: 

• X=l: duidt op het begin van het muziekstuk 

• T: de titel van het niuziekfragnient 

• M: de maatsoort 

• K=C: de toonaard is C: er staat dus geen jj of b vooraan 

• L=l/16: de basisnootduur is :j^ 

De functie 'schrijfabc' van TabuMusic genereert als output een .abc-file. Er zijn al- 
lerlei toepassingen beschikbaar die een .abc-files kunnen omzetten naar een .mid-file 
en tegelijk bladmuziek genereren in .pdf. Een voorbeeld van een gratis beschikbaar 
programma met deze functionaliteit is iABC. Dit opensource-programma is gratis te 
downloaden via http://abc.sourceforge.net/iabc/ (Walshaw, 2005), het staat 
tevens op de cd-rom. 

Nu we de working van TabuMusic besproken hebben, gaan we het in het volgende 
hoofdstuk verder testen. 



5 



Experimenten 



/ been searchin' . . . 
Searchin' . . . oh yeah 
Searchin' every which way 



(Leiber en StoUer, 1957) 



5.1 Inleiding 

In dit laatste hoofdstuk gaan we TabuMusic evalueren. We staan eerst stil bij de 
wegingscoefficienten van de doelfunctie. Vervolgens vergelijken we resultaten van 
melodische en ritniisclie optinialiseringen. Tenslotte volgt een persoonlijke niening 
van de gegenereerde niuziek. Alle testruns gebeurden met een Intel pentiuni IV 
met een kloksnellieid 1.4GIIz en 512MB RAM-gelieugen. Deze pc draait onder het 
operating system Windows XP. 

5.2 Wegingscoefficienten 

De doelfunctie is cruciaal bij het evalueren van muziek. Aangezien zij bestaat uit 
zes verschillende criteria (zie sectie 4.4), is het belangrijk om na te gaan hoe zwaar 
deze termen in de doelfunctie doorwegen. De doelfunctie is hieronder nog eens 
weergegeven: 



Icriterium 1 — 0. 27[ , Icriterium 2—0.011 , [criterium 3—0.491 , Icriterium 4 — 0.401 

score = aiJ ^j-^ l + ^^^ ^^ ^ + as^ 0:06 + ^'^ oTi 

Icriterium 5— 0.24| . |criteriuni 6 — 0.32| 
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/■ Criterium 1 ■/ 6,64 

/- — Criterium 2 — -/ 49,50 

A Criterium 3 ■/ 8,50 

/■ Criterium 4 ■/ 5,45 

/■ Criterium 5 ■/ 10,86 

/■ Criterium 6 ■/ 5,61 

Figuur 5.1: Mininiale en niaxiniale waarde van de 6 ternien in de doelfunctie. 

met flj = wegingscoefficient van criterium i 

In figuur 5.1 is de minimale en maximale waarde van de termen berekend (uiteraard 
zonder met de wegingscoefficienten rekening te houden). Deze waarden kunnen we 
berekenen aangezien criterium x steeds een waarde is tussen en 1. 

We merken dat bepaalde criteria overwegen. Om dit te compenseren kunnen we de 
wegingscoefficienten aanpassen. We doen dit op de volgende manier, waarbij a^ de 
wegingscoefficient van het z-de criterium is: 

49,50 ^7 ,r 



"■1 — 


6,64 


1 , '±U 


a2 


_ 49,50 
^ 49,50 


= 1 


0-3 = 


49,50 
8,50 


5,82 


a^ = 


49,50 
5,45 


9,08 


as = 


49,50 _ 
10,86 ^ 


4,56 


ag = 


49,50 _ 
5,61 


8,82 



In de rest van dit hoofdstuk gebruiken we steeds deze wegingscoefficienten. De cor- 
rectheid van deze coefficieiiten zou nog verder onderzocht kunnen worden. Misschien 
zijn bepaalde criteria belangrijker dan andere? Wij gaan er van uit dat ze allemaal 
even belangrijk zijn, maar dit zou nog verder onderzocht moeten worden. 
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5.3 Verbeteringen in score 

In deze sectie willen we evalueren hoe good het ontworpen tabu search- algoritme 
werkt. Daalt de score snel? Blijft ze dalen? We bekijken eerst afzonderhjk de 
nielodie en het ritnie. Daarna laten we het prograninia afwisselend ritnie en nielodie 
optimahseren. De tabu tenure bedroeg 1000 voor de 50 en 100 maten en 500 voor 
20 maten. 

5.3.1 Melodische verbetering 

In deze sectie vragen we ons of hoe groot de verbetering van de score is wanneer we 
enkel de nielodie optimahseren. Geraken we steeds aan een lage waarde? Vinden 
we snel een lage waarde? We hebben dit eerst onderzocht met een beperkt aantal 
iteraties (of moves). Nadien hebben we het programma wat langer laten draaien en 
de resultaten van 500 iteraties geanalyseerd. 

50 iteraties 

We hebben TabuMusic een aantal keer gerund. Telkens hebben we 50 moves ge- 
daan. Het gemiddelde voor van de testcases is voorgesteld in figuur 5.2. Er is een 
onderscheid gemaakt tussen muziekstukken van 20, 50 en 100 maten. 

Op de figuur is een duidelijk dalende trend te zien. De daling lijkt in het begin 
sterker te zijn dan op het einde. Verder valt ook op dat kortere muziekfragmenten 
gemiddeld een lagere beginscore hebben en sneller dalen. 

500 iteraties 

In figuur 5.3 zijn de gemiddelden voorgesteld van enkele test-fragmenten. We hebben 
TabuMusic nu 500 moves laten maken. Er is weer een onderscheid gemaakt tussen 
muziekstukken van 20, 50 en 100 maten. 

We zien dat er een sterke daling is in het begin. Na ongeveer 25 iteraties worden 
er geen grote dalingen meer gerealiseerd. Doordat we een zeer grote neighbourhood 
hebben kan het zijn dat we iets te lokaal blijven zoeken en daarom niet uit het lokale 
optimum geraken. Daarbij komt nog dat het ritme niet geoptimaliseerd wordt en de 
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Figuur 5.2: Resultaten van de melodische optimalisatie, 50 iteraties 
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Figuur 5.3: Resultaten van de melodische optimalisatie, 500 iteraties 
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score dus nooit kan worden. Oni uit het lokaal optiniuni te ontsnappen is er een 
diversificatie-mechanisme no dig. 

5.3.2 Ritmische verbetering 

In dit deel vragen we ons af hoe groot de verbetering in score is wanneer we enkel 
het ritme optinialiseren? Geraken we steeds aan een lage waarde? Geraken we er 
snel aan? We hebben net zoals bij de melodie een opsplitsing gemaakt tussen 50 
en 500 iteraties. De tabu tenure bedroeg 100 voor 50 en 100 maten en 10 voor 20 
maten. De lijst is dus minder lang als bij de melodie. Dit is zo gekozen omdat 
de tabulijst van het ritme iets stronger is, zie sectie 4.6. Ze zal dus sneller moves 
verbieden. Een voorbeeld maakt dit iets duidelijker: 

Op de tabulijst van de melodie kan bijvoorbeeld staan dat noot 2 met 
als midi-waarde 60 vervangen is. Daartegenover, op de tabulijst van het 
ritme zal enkel staan dat noot 2 is aangepast en niet dat dit bijvoorbeeld 
een achtstenoot was. In het geval van de melodie zal men noot 2 nog mogen 
aanpassen, zolang ze geen midi-waarde 60 krijgt. De ritmeoptimalisatie is 
iets strenger en zegt dat er niets meer mag veranderen aan noot 2 zolang 
ze op de tabulijst staat. 

Door een kortere tabulijst te nemen voorkomen we dat er na een tijdje geen enkele 
move meer mogelijk is. 

50 iteraties 

In figuur 5.4 zijn de resultaten voorgesteld, analoog aan de vorige figuren. Nu is 
echter het ritme geoptimaliseerd. 

Het valt op dat er een sterk dalende trend is bij korte muziekfragmenten. Bij langere 
stukken is dit echter niet het geval. Het zou kunnen dat er meer dan 50 moves nodig 
zijn om het ritme van een lang muziekstuk te herschikken. Het zou ook kunnen dat 
er een zeer snelle convergentie is. Het is daarom nuttig om het programma iets meer 
moves te laten uitvoeren. 
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Figuur 5.4: Resultaten van de ritmische optinialisatie, 50 iteraties 

500 iteraties 

In figuur 5.5 zijn de resultaten voorgesteld, analoog aan de vorige figuren. Er is nu 
geoptimaliseerd met 500 iteraties 

Er is een continue daling van de doelfunctie. De langere niuziekstukken convergeren 
dus niet zeer snel, ze hebben iets nieer tijd nodig oni te zakken. Het valt ook op dat 
er al iets nieer rekentijd voor nodig is (meer dan 5 niinuten). Kortere fragnienten 
dalen wel snel en blijven dan nog traag verder dalen. Wanneer we het progranima 
nog veel langer laten lopen komt er misschien wel een punt van convergentie, dit 
zou nog verder getest kunnen worden. We willen echter niet alleen ritmisch goede 
muziek, ook de nielodie is belangrijk. Daaroni gaan we hierna analyseren wat er 
gebeurd als we beiden optimaliseren. 

5.3.3 Ritme en melodie afwisselend verbeterd 



Dit deel is eigenlijk het belangrijkste. We vragen ons af hoe groot de verbetering 
van de score is wanneer we zowel het ritme als de melodie optimaliseren? Geraken 
we steeds aan een lage waarde? Geraken we er snel aan? We hebben steeds 10 
moves het ritme geoptimaliseerd, en dan 10 moves de melodie, enzoverder. 
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Figuur 5.5: Resultaten van de ritmische optimalisatie, 500 iteraties 



50 iteraties 

In figuur 5.6 zijn de gemiddelde waarden uitgezet tegen het aantal iteraties. We 
niaken weeroni een opsplitsing tussen 20, 50 en 100 maten. 

AUe curves dalen. Dit was te verwachten aangezien de melodisclie en ritmische 
optimalisatie elk apart ook daalden in het begin. Ook hier zien we, net zoals bij het 
ritme, dat korte fragmenten sneller dalen. 



500 iteraties 

In figuur 5.7 zijn de gemiddelde waarden uitgezet tegen het aantal iteraties. We 
maken weerom een opsplitsing tussen 20, 50 en 100 maten. 

Ook na 500 moves merken we dat de score nog steeds blijft dalen. We herinneren 
ons dat dit niet het geval was wanneer we enkel de melodie verbeterden. Voor grote 
fragmenten van 100 maten is de daling zeer zwak. Het zou kunnen dat ze nog meer 
tijd no dig heeft om sterk te dalen omdat de te verkennen ruimte veel groter is. 
Eerder zagen we bij het ritme, dat grotere fragmenten meer tijd nodig hadden om 
een betere score te krijgen. De rekentijd is daarenboven ook langer bij 100 maten 
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Figuur 5.6: Resultaten van de geconibineerde optimalisatie, 50 iteraties 
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De gemiddelde rekentijd bedroeg: 



400 500 

10 sec voor 20 maten 
46 sec voor 50 maten 
197 sec voor 100 maten 



Figuur 5.7: Resultaten van de geconibineerde optimalisatie, 500 iteraties 
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(ongeveer 3.5 min). Het is nuttig om enkele tests te doen met zeer veel iteraties om 
te zien of de score van grotere stukken nog blijft dalen. 

Nog meer iteraties 

Door de lange rekentijd die nodig is om veel iteraties te doen, hebben we de testruns 
zelf voorgesteld in figuur 5.8 in plaats van het gemiddelde te nemen van meerdere 
tests. 



CD 
O 

o 




n 


Test A - 


-50 


maten 


♦ 


TestB- 


-50 


maten 


V 


TestC 


-50 


maten 



200 400 800 1200 1600 2400 

Aantal iteraties 



De gemiddelde rekentijd bedroeg: 5 min voor 50 maten 

Figuur 5.8: Resultaten van de gecombineerde optimalisatie, 2400 iteraties 



We merken dat er vooral in het begin een zeer sterk dalende trend is. Nadien is 
de daling lets minder sterk, test B blijft ongeveer constant. Het kan zijn dat ze 
na een nog langere tijd verder dalen. Het is echter meer waarschijnlijk dat ze na 
2400 moves ongeveer gaan convergeren. Dit kan te wijten zijn aan de aard van 
de gekozen moves. De neighbourhood is telkens zeer groot en de inividuen uit de 
neighbourhood verschillen zeer weinig van de originele muziek. Het kan daarom 
interessant om bijvoorbeeld als move te nemen: wissel twee noten om in plaats van 
een noot te veranderen. 
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5.4 Algemene muzikale kwaliteit 

Een voorbeeld van een gegenereerd muziekstuk is gegeven in bijlage E. Er bestaat 
geen twijfel dat de willekeurig gegenereerde niuziek edit willekeurig is. Deze muziek 
is niet aangenaam oni naar te luisteren en klinkt zeer 'random'. Het is de persoonlijke 
niening van de auteur dat geoptinialiseerde niuziek, met een lage score, iets beter 
klinkt dan de willekeurige. Er zitten minder grote intervals in. De overgangen 
klinken vaak beter en liet ritnie is consistenter. 

De gegenereerde muziek is ecliter maar zo goed als de doelfunctie toelaat. Er is geen 
enkele criterium dat iets zegt over de contour, een tlienia, of dergelijk. Tocli lijkt 
de geoptinialiseerde muziek reeds goed te klinken, ondanks de beperkingen van de 
doelfunctie. 



5.5 Conclusie 

Tabu search kan zeker ook worden aangewend op het vlak van computer assisted 
composing. TabuMusic slaagt er redelijk goed in om de score te minimaliseren. 
We vinden snel een redelijk lage score. Na een tijdje zijn de dalingen ecliter niet 
meer zo groot begint het algoritme te convergeren. Dit kan liggen aan het feit dat 
de moves zeer beperkt zijn. In verder onderzoek zou men meer en andere moves 
kunnen inbouwen, om zo meer differentiatie te krijgen. 

De muziek die gegenereerd wordt klinkt niet slecht, rekening houdend met de be- 
perkingen van de doelfunctie. Het zou zeer interessant zijn om de doelfunctie verder 
uit te breiden zodat ze ook toelaat dat er een thenia inkonit. We zouden ook de 
optimale waarden van de verschillende criteria na kunnen gaan. Dit kunnen we doen 
door de waarde van elk criteria van bestaande bekende fragmenten te berekenen en 
hiervan het gemiddelde te nenien. 

Een andere interessante uitbreiding is het optimaliseren van polyfone muziek. In de 
hedendaagse klassieke muziek is het zeldzaam dat er maar een noot tegelijk klinkt. 
Ook zou het interessant zijn om bestaande muziekstukken in te geven en ze te 
optimaliseren. Zouden ze nog verbeteren? 



Besluit 



Muziek en computers gaan zeer goed samen. Sinds liet begin van de 20ste eeuw 
is er veel gebeurd op liet gebied van computer assisted composing of CAC. Steeds 
krachtigere computers laten toe om ingewikkeldere programma's te schrijven die 
beter met muziek om kunnen. 

We hebben gezien dat we het componeren van muziek kunnen beschouwen als een 
optimaliseringsprobleem. Daardoor kunnen we teclinieken uit de wereld van de op- 
timalisatie gebruiken om het compositieproces te versnellen. Metaheuristieken zijn 
intelligente zoekstrategieen die hierbij van pas kunnen komen. Genetische algorit- 
men, een bepaalde soort metaheuristiek, is een veelgebruikte techniek bij het ont- 
werpen van CAC-programma's. De besproken projecten GenJam en SIEDP, werken 
allebei op basis van een genetisch algoritme. Zij zijn voorbeelden van een succesvolle 
toepassing van metaheuristieken voor CAC. 

Van een andere metaheuristiek, namelijk tabu search, zijn er weinig of geen voorbeel- 
den bekend. Het leek ons daarom interessant om zelf een toepassing te ontwikkelen 
(TabuMusic). Het ontworpen programma is bijgeleverd op de cd-rom en is in staat 
om monofone muziekfragmenten te componeren. Enkele voorbeelden van muziek 
gegenereerd door TabuMusic staan eveneens op de cd-rom. TabuMusic levert reeds 
goede resultaten, toch is er na een langere tijd convergentie. Een uitbreiding, die 
enkele mechanismen voor differentie voorziet, zou hiervoor een oplossing kunnen 
bieden. 

Tabu Search is een zoekstrategie die meer gerichter is dan de veelgebruikte gene- 
tische algoritmen. Deze local search strategie komt normalerwijze niet in cirkels 
vast te zitten door het gebruik van geheugenstructuren zoals de tabulijst. Voor 
het verkennen van de enorme ruimte van mogelijke muziekstukken is dit zeker han- 
dig. We kunnen concluderen dat tabu search een nuttig hulpmiddel kan zijn bij het 
ontwerpen van CAC-toepassingen. 
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DEEL III 



Bijlagen 



A 



Textuele weergave van een 
MIDI-bestand 



Dit voorbeeld is gemaakt met behulp van (Glatt, 1998). 
MThd I Format=l I # of Tracks=2 I Division=96 

Time Event 
1: 1: I Tempo I BPM=60 I micros\quarter=1000000 
I Time Sig I 4/4 I MIDI-clocks\click=24 I 32nds\quarter= 
I Key Sig I C Major I 
I End of track I 

Time Event 
1:1: I On Note I chan= 1 I pitch=E 4 I vol=80 

3: I Off Note I chan= 1 I pitch=e 4 I vol=64 

I (Off) Note I chan= 1 I pitch=c 3 
24 I Off Note I chan= 1 I pitch=c 3 I vol=64 
I On Note I chan= 1 I pitch=G 3 I vol=80 

4: I Off Note I chan= 1 I pitch=g 3 I vol=64 
I On Note I chan= 1 I pitch=G 3 I vol=80 

2: 1: 72 I Off Note I chan= 1 I pitch=g 3 I vol=64 
I On Note I chan= 1 I pitch=F#4 I vol=80 

4: 24 I Off Note I chan= 1 I pitch=f#4 I vol=64 
I On Note I chan= 1 I pitch=G 4 I vol=80 

3: 2: 72 I Off Note I chan= 1 I pitch=g 4 I vol=64 
I On Note I chan= 1 I pitch=D 3 I vol=80 
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4: 48 I Off Note 

I On Note I chan= 1 
4: 4: I Off Note 

I On Note I chan= 1 
5: 2: 48 I Off Note 

I On Note I chan= 1 
6:2: I Off Note 

I On Note I chan= 1 

I End of track I 



chan= 1 I pitch=d 3 I vol=64 

I pitch=E 3 I vol=80 
chan= 1 I pitch=e 3 I vol=64 

I pitch=B 4 I vol=80 
chan= 1 I pitch=b 4 I vol=64 

I pitch=F#3 I vol=80 
chan= 1 I pitch=f#3 I vol=64 

I pitch=C 4 I vol=80 



B 



Een Orchestra- en een 
Score-bestand van CSound 



Voorbeeld overgenomen van (FEMS/CAMIL, 2004). 



Orchestra-bestand 



; basicgrain . ore 





70_01_l.orc 


; synthesis: 


granular 




basic granular synthesis 




uniform pseudorandom selection of 




grain frequency, inter-grain time, 




and grain length 


; coded: 


6/95 RKP 


sr=48000 




kr=480 




ksmps=100 




nchiils=l 




instr 1 




lamp = p4 


; amplitude of overall event 


ibpch = p5 


;basepch (center frequency) 
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ibw = p6 
of grains. 



imiigt 


= P7 


imxigt 


= p8 


imigl 


= p9 


imxgl 


= plO 


iseed 


= pll 


igefn 


= pl2 



; bandwidth of random frequency variation 

;1 = one semitone variation above and 

below the base frequency 
; minimum inter-grain time 
; maximum inter-grain time 
; minimum grain length 
; maximum grain length 
;seed value for randh 
; grain envelope function number 



; inter-grain time: uniform pseudorandom selection 
between imiigt and imxgt 

krl randh (imxigt-imiigt)/2, 1/imiigt , iseed 
kigt = imiigt+ (krl+( (imxigt-imiigt) /2) ) 

; grain length: uniform pseudorandom selection 
between imigl and imxgl 

kr2 randh (imxgl-imigl) /2, 1/imiigt , iseed 
kgl = imigl+(kr2+((imxgl-imigl)/2)) 

; grain cps random variation: uniform pseudorandom 

selection. 

;mean = ibpch, range = (ibpch) + - ibw 

kvar randh ibw/12, 1/imiigt , iseed 

; envelope for overall event 

kenv linen iamp,p3* . 1 ,p3,p3* .3 



reset : 



contin: 



timout 0, i (kgl)+i (kigt) , contin 

reinit reset 

andx line 0, i (kgl) ,4096 

agenv tablei andx, igefn 

ioct = octpch(ibpch) 

icps = cpsoct (ioct+i(kvar) ) 

al oscili agenv, icps, 1 
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endin 



Score-bestand 



out al*i(kenv) 



;basicgraiii. sco 

;70_01_l.sco 

; coded: RKP 6/95 

;iamp amplitude of overall event 
;ibpch base pch (center frequency) 
; ibw bandwidth of random frequency 
variation of grains, 

;1 = one semitone variation above 
and below the base frequency 
;imiigt minimum inter-grain time 
; imxigt maximum inter-grain time 
; imigl minimum grain length 
; imxgl maximum grain length 
;iseed seed value for randh 
; igef n grain envelope function number 



; GEN functions 
; waveform 
fl 4096 
grain production 



10 



;sine function for 



; grain envelope 

f31 4096 7 



1024 1 2048 1 1024 



; score 



; ins strt dur lamp ibpch ibw imiigt 



no 



imxigt 


imigl 


imxgl 


iseed 


igefn 






il 





1 


10000 


9.09 


1 


.02 


.1 


.02 


.05 


.001 


31 






s 
; ins 


strt 


dur 


lamp 


ibpch 


ibw 


imiigt 


imxigt 


imigl 


imxgl 


iseed 


igefn 






il 





1 


10000 


9.09 


2 


.02 


.1 


.02 


.05 


.001 


31 






s 
; ins 


strt 


dur 


lamp 


ibpch 


ibw 


imiigt 


imxigt 


imigl 


imxgl 


iseed 


igefn 






il 





1 


10000 


9.09 


3 


.02 


.1 


.02 


.05 


.001 


31 






s 
; ins 


strt 


dur 


lamp 


ibpch 


ibw 


imiigt 


imxigt 


imigl 


imxgl 


iseed 


igefn 






il 





1 


10000 


9.09 


6 


.02 


.1 


.02 


.05 


.001 
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Broncode TabuMusic 



uses crt , sysutils; 

const 

maxtijden = 1600; 

maxmaten = 100; 

aantaltonen = 14; {aantaltonen in een toonladder, 

wanneer we 2 octaven nemen van do to si} 
aantalnoten = 24; 
laagstenoot = 60; 
maxtabu = 1000; 
maxtabu2 = 1000; 



type 

muztype = array[1..2, 1 . .maxtijden] of longint; 
toontype = array [1 .. aantaltonen] of longint; 
tabutype = array [1.. 2, 1.. maxtabu] of longint; 
tabu2type = array [1 . .maxtabu2] of longint; 

var 

uur, min, sec, ms : word; 

infile, logfile : text; 

muziek, muziektijdelijk2 : muztype; 

toonladder : toontype; 

tabulijst : tabutype; 

tabulijst2 : tabu2type; 
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aantaltabu, aaiitaltabu2 : longint; 

iters, i, j, k, kies, show, somi j , somik, a, b, d, e, f, z, s, y, noot, 

consec, positie, aantaltijden, grondtoon, aantalmaten : longint; 
hulp, punten : real; 

wegingl , weging2, wegingS, weging4, wegingS, weging6 : real; 
bestandsnaam, bestandsnaamin, bestandsnaammuzl , bestandsnaammuz2, 

bestandsnaamlog : string; 

start, einde, verschil, totaal : TDateTime; 

error : boolean; 

aspiratiescore : real; 

aspiratie : muztype; 

function doelfunctie(muziekd : muztype) : real; 

var 

crit5, critl, critS, critlO, scorel, crit7, crit2, diss, critS, 

critll, crit4, crit6 : real; 
aantal: longint; {aantal verschillende tonen} 
max, min, herhaald, c, aantalntoon, somb,geenrust , aantalrust, hoogste, 

laagste, duren, ok, b, aantalzrust : longint; 
intervals, stijgend, verschil : longint; 

{stijgend = aantal stijgende intervallen} 

noten: longint ; {aantal noten waaruit het stuk bestaat} 

notenlijst : array [1 .. aantalnoten] of longint; 

muzrust , muziekr, muziek2 : muztype; 

duur : array [1.. 16] of longint; 

doell, doel2, doel3, doel4, doel5, doel6 : real; 

begin 

{eerste criterium: pitch variety= aantal verschillende noten 

(=aantal)/aantalnoten(=noten)} 
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{1. de verschillende noten tellen} 
{lijst leegmaken} 
for i:=l to aantalnoten do 
notenlijst [i] := 0; 

{notenlijst = 1 als de noot erin zit} 
for i:=l to aantaltijden do 
begin 

k := (muziekd[l , i] -laagstenoot+1) ; 
if (k <> (-laagstenoot)) 

then notenlijst [k] := 1; {rusten(-l) niet meetellen} 
end; 



{aantal enen tellen in notenlijst} 

aantal := 0; 

for i:=l to aantalnoten do 

begin 

if notenlijst [i] = 1 

then aantal := aantal +1; 

end; 



{Uit hoeveel noten bestaat het stuk?} 
noten := 0; 

for i:=l to aantaltijden do 
begin 

if (muziekd[2,i] = 0) and (muziekd [1 , i] <>(-!) ) 
then noten := noten +1; 
end; 



critl :=(aantal/noten) ; 



{criterium 2: dissonant intervals. Een interval van 10 krijgt een 
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dissonantiescore (=diss) van 0.5. Intervallen van 6, 11 of >= 13 

krijgen een diss van 1.0. 

De dissonantiescore wordt gedeeld door het aantal intervals} 

{allenoten in een matrix zetten zonder de rusten en lange noten: muzrust} 
aantalzrust := 0; 
for i:=l to aantaltijden do 
begin 
if (muziekd[l ,i] <> -1) and (inuziekd[2, i] = 0) then 
begin 

aantalzrust := aantalzrust+1 ; 
muzrust [1 , aantalzrust] := muziekd [1 , i] ; 
muzrust [2, aantalzrust] := muziekd [2 , i] ; 
end; 
end; 



{bepalen dissonantiescore} 

diss :=0; 

for i:=l to (aantalzrust-1) do 
begin 
verschil := ABS (muzrust [1 , i] -muzrust [1 , i+1] ) ; 



= diss + 0.5; end; 
= diss +1; end; 
= diss +1; end; 



if (verschil=10) then begin diss 
if (verschil=6) then begin diss 
if (verschil=ll) then begin diss 
if (verschil>=13) then begin diss := diss +1; end; 
end; 
if noten = 1 then crit2 := 1 

else crit2 := diss/(noten-l) ; 



{criterium 3: contour direction. De som van alle stijgende intervals 
gedeeld door de som van alle intervals} 
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intervals :=0; 
stijgend := 0; 

for i:=l to (aantalzrust-1) do 
begin 
if muzrust [1 , i+1] > muzrust[l,i] 

then stijgend := stijgend + muzrust [1 , i+1] -muzrust [1 , i] ; 
intervals := intervals + ABS (muzrust [1 , i+1] - muzrust [1 , i] ) ; 
end; 
if intervals = then crit3 := 1 

else critS := stijgend/ (intervals) ; 



{criterium 4: contour stability, aantal opeenvolgende 

intervallen die in dezelfde richting evolueren. } 

{tellen aantal opeenvolgende intervals in dezelfde richting} 

consec := 0; 

for i:=l to (aantalzrust-2) do 

begin 

if (muzrust [l,i+2] >= muzrust [1 , i+1] ) and 

(muzrust [1 , i+1] >= muzrust [1 , i] ) 
then consec := consec + 1; 
if (muzrust [l,i+2] <= muzrust [1 , i+1] ) and 

(muzrust [1 , i+1] <= muzrust [1 , i] ) 
then consec := consec + 1; 
end; 
if noten=2 then crit4 := 1 
else crit4 := consec/ (noten-2) ; 



{criterium 5: rythmic variety. 

aantal verschillende nootduraties gebruikt / 16 

(aantal verschillende nootduraties ongeveer mogelijk} 
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for i:=l to 16 do 
duur[i] := 0; 

i:=l; 

while i<= aantaltijden do 

begin 
b := 0; 

if muziekd [2, i] = 1 then 
begin 

while (muziekd [2, i+b] = 1) do 

begin 
{gewoon tellen} 
b := b +1; 
end; 

end; 

if muziekd [2, i+b] = then 
begin 

somb := b+1; 

duur[somb] := duur [somb] +1 ; 
end; 

i := i+b+1; 



end; 



{aantal verschillende tijdsduren tellen} 
duren := 0; 
for i := 1 to 16 do 
begin 
if duur[i] <> 

then duren := duren + 1; 
end; 
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crit5 := duren/16; 



{Criterium 6: rythmic change 
(Max noot duratie + min noot duratie) / 16} 
{zoek max noot duratie: max} 
max : = 1 ; 
min := 1; 
ok := 0; 
i:=16; 

while ok = do 
begin 

if duur[i] <> then 
begin 
max : = i ; 
ok := 1; 
end; 
i := i-1; 
end; 

{vind de minimum duur} 
ok := 0; 
i := 1; 

while ok = do 
begin 

if duur[i] <> then 
begin 
min := i; 
ok := 1; 
end; 
i := i+1 
end; 



crit6 := (max-min) /16 ; 
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doell 
doel2 
doel3 
doel4 
doel5 
doel6 



= (ABS(critl-0.27))/0.11 
= (ABS(crit2-0.01))/0.02 
= (ABS(crit3-0.49))/0.06 
= (ABS(crit4-0.40))/0.11 
= (ABS(crit5-0.24))/0.07 
= (ABS(crit6-0.32))/0.11 



{totale score} 

doelfunctie := (wegingl*doell)+(weging2*doel2)+(weging3*doel3)+ 

(weging4*doel4)+(weging5*doel5)+(weging6*doel6) ; 
if show = 1 then 

begin 

writelnClogfile, 'Pitch variety: ', #9, #9, doell:0:6); 

writelndogf ile, 'Dissonant intervals:', #9, doel2:0:6); 

writelnClogfile, 'Contour direction: ', #9, doel3:0:6); 

writelnClogfile, 'Contour stability: ', #9, doel4:0:6); 

writelnClogfile, 'Rhythmic variety: ', #9, doel5:0:6); 

writelnClogfile, 'Rhythmic range: ', #9, doel6:0:6); 

end; 



end; 



procedure schrijf abcCmuziekschrijf : muztype ; 

bestandsnaamschrijf : string) ; 

var 

outfile : text; 

teller, tijd : longint; 

abc : string; 



begin 

assignCoutf ile, bestandsnaamschrijf) ; 
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rewrite(outf ile) ; 

writeln(outfile, 'X:l') ; 

writeln(outf ile, 'TiMuziekfragment TabuMusic') ; 

writelnCoutfile, 'M:4/4') ; 

writelnCoutfile, 'K:C') ; 

writelnCoutfile, 'L:l/16') ; 



tijd := 0; 

teller := 1; 

for i:=l to aantaltijden do 

begin 

{als 2e lijn = dan schrijven} 

if muziekschrijf [2,i] = then 

begin 

{noot schrijven naar bestand} 

if muziekschrijf [1 ,i] = 60 then write(outf ile 



if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 



if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 

if muziekschrijf [1 



i] = 61 then write(outf ile 

i] = 62 then write(outf ile 

i] = 63 then write(outf ile 

i] = 64 then write(outf ile 

i] = 65 then write(outf ile 

i] = 66 then write(outf ile 

i] = 67 then write(outf ile 

i] = 68 then write(outf ile 

i] = 69 then write(outf ile 

i] = 70 then write(outf ile 

i] = 71 then write(outf ile 



if muziekschrijf [1 ,i] = 72 then write(outf ile, 



i] = 73 then write(outf ile 

i] = 74 then write(outf ile 

i] = 75 then write(outf ile 

i] = 76 then write(outf ile 

i] = 77 then write(outf ile 

i] = 78 then write(outf ile 

i] = 79 then write(outf ile 

i] = 80 then write(outf ile 
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if muziekschrijf [1 ,i] = 81 then write(outfile, ' a' .teller, ' '); 

if muziekschrijf [1 ,i] = 82 then write(outfile, ' _b' , teller , ' '); 

if muziekschrijf [1 , i] = 83 then write(outfile, 'b' , teller, ' '); 

if muziekschrijf [1 ,i] = -1 then write(outfile, 'z' , teller, ' '); 

teller := 1; {we beginnen terug met een nieuwe noot die 

min 1 tel duurt} 

tijd := tijd + 1; 

if tijd = 12 then 

begin 

writeln(outf ile, ' ') ; writeln(outf ile, ' ') ; 

tijd := 0; 

end; 
end; 

if muziekschrijf [2, i] = 1 then 
begin 

teller := teller + 1 ; 
end; 

end; 

close(outf ile) ; 



end; 



function noottransponeren(muziektransp : muztype; 
scoretransp: real): muztype; 

var 

muziekmeer, muziektijdelijk : muztype; 
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scoretijdelijk, scoremeer, beginscore : real; 

okl , X , V, dalen, gelijk, tabuxmeer, tabuymeer: longint; 

vervangen : boolean; 

tabux, tabuy : longint; 

begin 

tabux := 0; 

tabuy := 0; 

noottransponeren := muziektransp; 

muziektijdelijk := muziektransp; 

dalen := 1; 

scoremeer := 10000; 

beginscore := scoretransp; 

for X := 1 to aantaltonen do 

begin 

{probeer elke noot te vervangen} 

y := 1; 

while y <= aantaltijden do 
begin 

{vervang 1 noot} 
z := 0; 
okl := 0; 
while okl = do 
begin 
muziektijdelijk [1 ,y+z] := toonladder [x] ; 

if muziektijdelijk[2,y+z] = then 
begin 
okl := 1; 
end; 
z := z+1; 
end; 
scoretijdelijk := doelf unctie(muziekti jdelijk) ; 
vervangen := false; 



123 



{geen betere gevonden} 
{best mogelijke slechte oplossing onthouden} 
if scoretijdelijk > beginscore then 
begin 

if scoremeer > scoretijdelijk then 
begin 

scoremeer := scoretijdelijk; 
muziekmeer := muziektijdelijk; 

tabuxmeer := y; 

tabuymeer := muziektijdelijk [1 ,y] ; 
end; 
end; 

{betere gevonden!} 

{is het dezelfde oplossing?} 

gelijk := 1; 
for V := 1 to aantaltijden do 
begin 
if muziektijdelijk[l ,v] <> muziektransp [1 ,v] then 

gelijk := 0; 
end; 
for V := 1 to aantaltijden do 
begin 
if muziektijdelijk [2, v] <> muziektransp [2, v] then 

gelijk := 0; 
end; 



if (scoretijdelijk <= scoretransp) and (gelijk = 0) then 
begin 
vervangen : = true ; 
for v:=l to aantaltabu do 
begin 

if (muziektijdelijk[l ,y] =tabulijst [2,v] ) then 
begin 
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if (tabulijst [l,v]=y) then 
begin 
vervangen:= false; 
end; 
end; 
end; 

{toch vervangen als het aspiratieniveau is bereikt} 

if scoretijdelijk < aspiratiescore then 

begin 

aspiratie := muziektijdelijk; 

aspiratiescore := scoretijdelijk; 

vervangen : = true ; 

end; 

{geen TABU} 

if vervangen = true then 
begin 
{noot vervangen} 
dalen := 0; 

scoretransp := scoretijdelijk; 
noottransponeren := muziektijdelijk; 
{onthouden welke getallen er in de tabulijst moeten komen} 
tabux := y; 

tabuy := muziektijdelijk [1 ,y] ; 
end; 
end; 

muziektijdelijk := muziektransp; 

y := y+z; 

end; 
end; 

if dalen = 1 then 
begin 
scoretransp := scoremeer; 
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noottransponeren := muziekmeer ; 
tabux := tabuxmeer; 
tabuy := tabuymeer; 
end: 



{tabulijst bijvullen als het is aangenomen} 
if (dalen = 0) or (dalen=l) then 
begin 

{tabulijst opschuiven} 
for V := (aantaltabu-1) downto 1 do 
begin 
tabulijst [1 ,v+l] := tabulijst [1 ,v] ; 
tabulijst [2, v+1] := tabulijst [2, v] ; 
end; 

{tabulijst bijvullen} 
tabulijst [1 , 1] := tabux; 
tabulijst [2, 1] := tabuy; 

end; 
end; 



function ritmeverbeterenCmuziekrverb : muztype; punten2:real) : muztype; 

var 

ritmeverb, muziekmeer2 : muztype; 

telnoten: longint; 

cl, rest, x2, yl , v2, u, dalen2, gelijk2: longint; 

ritmemuziek, ritmeinuziek2 : array [1.. 2, 1 . .maxtijden+1] of longint; 

vervangen2, ok2 : boolean; 

scoretijdelijk2, scoremeer2, beginscore2 : real; 

tabuxl, tabuxlmeer : longint; 



begin 
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ritmeverb := muziekrverb; 
muziektijdelijk2 := muziekrverb; 
ritmeverbeteren := muziekrverb; 
dalen2 := 1; 
scoremeer2 := 10000; 
beginscore2 := punten2; 
tabuxl := 0; 



{noten tellen} 

telnoten := 0; 

for x2 := 1 to aantaltijden do 

begin 

if muziekrverb [2, x2] = then 

begin 

telnoten := telnoten + 1; 

end; 
end; 



{ritmemuziek initialiseren} 

for x2 := 1 to 2 do 

begin 

for yl := 1 to (maxtijden+1) do 

ritmemuziek [x2,yl] := 0; 

end; 



{in speciale matrix zetten. Rij : midiwaarde en Kolom 
aantaltijden dat de noot duurt} 

v2 := 0; 

u := 1; 

for x2 := 1 to aantaltijden do 

begin 

if muziekrverb [2, x2] = 1 then 
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begin 

v2 := v2 +1; 
end; 

if muziekrverb [2,x2] = then 
begin 

v2 := v2 + 1; 

ritmemuziek [1 ,u] := muziekrverb [1 ,x2] ; 

ritmemuziek [2,u] : = v2 ; 

u := u + 1; 
v2 := 0; 
end; 



end; 



{probeer het ritme van de noot te vervangen door elk mogelijk ander} 
{eerst voor de vde noot noot} 
for v2 := 1 to telnoten do 
begin 

for u := 1 to 16 do 
begin 

ritinemuziek2 := ritmemuziek; 

ritmemuziek2 [2, v2] := u; 

{terug omzetten in goeie matrix} 

yl := 1; 
cl := 1; 
rest := ritmemuziek2 [2, cl] ; {rest: aantal tijden 

dat de noot nog verderklinkt} 
while yl <= aantaltijden do 
begin 
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ok2 := true; 

{per kolom van ritmemuziek: } 

{begin} 

muziektijdelijk2 [1 ,yl] := ritmemuziek2 [1 , cl] ; 
{als er geen noten meer zijn: } 
if ritmemuziek2 [2, cl] = then 
begin 

muziekti jdelijk2 [1 ,yl] := ritmemuziek [1 ,v2] ; 
muziektijdelijk2 [2,yl] := 1; 
rest := 0; 
end; 

if rest = 1 then 
begin 

muziektijdelijk2 [2,yl] := 0; 
cl := cl + 1; 

rest := ritmemuziek2 [2, cl] ; 
ok2 := false; 
end; 

if (rest > 1) and (ok2) then 
begin 

muziektijdelijk2 [2,yl] := 1; 
rest := rest - 1; 
end; 

yl := yl + 1; 
end; 



{laatste noot moet stoppen: } 
muziekti j deli jk2 [2, aantaltij den] := 0; 
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{doelfunctie berekenen} 

scoretijdelijk2 := doelfunctie (muziektij deli jk2) ; 
vervangen2 := false; 

{geen betere gevonden} 

{best mogelijke oplossing onthouden} 

if scoretijdelijk2 > beginscore2 then 

begin 

if scoremeer2 > scoretijdelijk2 then 

begin 

scoremeer2 := scoretijdelijk2; 
muziekmeer2 := muziektijdelijk2; 

tabuxlmeer := v2; 
end; 
end; 

{is het dezelfde oplossing?} 
gelijk2 := 1; 
for x2 := 1 to aantaltijden do 
begin 
if muziektijdelijk2 [1 ,x2] <> muziekrverb [1 ,x2] then 
gelijk2 := 0; 
end; 
for x2 := 1 to aantaltijden do 
begin 
if muziektijdelijk2 [2,x2] <> muziekrverb [2, x2] then 

gelijk2 := 0; 
end; 



{lage score} 
if (scoretijdelijk2 <= punten2) and (gelijk2 = 0) then 
begin 

vervangen2 := true; 
for x2 := 1 to aantaltabu2 do 
begin 
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if (tabulijst2[x2]=v2) then 
begin 

vervangen2 := false; 
end; 
end; 
end; 

{toch vervangen als het aspiratieniveau is bereikt} 

if scoretijdelijk2 < aspiratiescore then 

begin 

aspiratie := muziektijdelijk2; 

aspiratiescore := scoretijdelijk2; 

vervangen2 : = true ; 

end; 

{hogere score} 
if vervangen2 = true then 

begin 

{nieuwe oplossing aannemen} 
dalen2 := 0; 

punten2 := scoretijdelijk2; 

ritmeverbeteren := muziektijdelijk2; 
{onthouden getallen tabulijst} 
tabuxl := v2 ; 

end; 
end; 
end; 

if dalen2 = 1 then 

begin 

punten2 := scoremeer2; 

ritmeverbeteren := muziekmeer2; 

tabuxl := tabuxlmeer; 
end; 

if (dalen2 = 0) or (dalen2=l) then 



131 



begin 
{tabulijst} 

for x2 := (aaiitaltabu2-l) downto 1 do 
begin 

{tabulijst opschuiven} 
tabulijst2[x2+l] := tabulijst2 [x2] ; 
end; 
{tabulijst bijvullen} 
tabulijst2 [1] := tabuxl; 
{tabulijst drukken} 
end; 

end; 



begin 

{scherminstellingen} 
clrscr ; 

textbackground (white) ; 
textcolor(black) ; 



{tabulijst initialiseren} 
for i := 1 to 2 do 

begin 

for j := 1 to aantaltabu do 
tabulijst [i,j] := 0; 

end; 

for i := 1 to aantaltabu2 do 
begin 

tabulijst2[i] := 0; 
end; 

{initialisatie} 
error := true; 
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{init infile} 

bestandsnaamin := 'inputmuziek.txt'; 
assign (infile, bestandsnaamin); 
reset (infile) 



readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 
readln(inf ile 



grondtoon) ; 
aantalmaten) ; 
bestandsnaammuzl) ; 
bestandsnaaminuz2) ; 
bestandsnaamlog) ; 
iters) ; 
aantaltabu) ; 
aantaltabu2) ; 
wegingl) 
weging2) 
wegingS) 
weging4) 
wegingS) 
weging6) 



{init logfile} 

assigndogf ile, bestandsnaamlog) ; 

rewritedogf ile) ; 

totaal := 0; 



if (grondtoon <60) or (grondtoon>71) then 
begin 

writelnCU hebt een verkeerde grondtoon opgegeven! ' ) ; 
error := false; 
end; 



if (aantalmaten > maxmaten) then 
begin 
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writelnCU hebt een te groot aantal maten opgegeven ! ' ) ; 
error := false; 
end; 



{vul toonladderverzameling} 

{zoek de laagste noot van de toonladder} 

noot := grondtoon-12; 

positie := 1; 

while noot < laagstenoot do 

begin 

if (positie = 1) or (positie = 2) or (positie = 4) 

or (positie = 5) or (positie = 6) 

then noot := noot + 2 

else noot := noot + 1; 
positie := positie +1; 
if positie = 8 then positie := 1; 
end; 
toonladder [1] := noot; 

{genereer de rest van de toonladder} 
for i:=2 to aantaltonen do 

begin 

if (positie = 1) or (positie = 2) or (positie = 4) 

or (positie = 5) or (positie = 6) 

then noot := noot + 2 

else noot := noot + 1; 
positie := positie +1; 
if positie = 8 then positie := 1; 
toonladder [i] := noot; 
end; 



{genereer een eerste willekeurige oplossing} 
aantaltijden := aantalmaten * 16; 
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randomize ; 

i := 1; 

while i<= aantaltijden do 

begin 

a := random (aant alt onen+1 ) ; {random toon} 

k := random (1 6) +1 ; {random duur} 

if ((i+k) > aantaltijden) then 
begin 
k := aantaltijden - i+l;{-k+l} 
end; 



{eerste rij opvullen} 
for j := 1 to k do 
begin 

if a=0 then 
begin 

somij := i+j-1; 
muziek[l jSomij] := -1; 
end; 

if a<> then 

begin 

somij := i+j-1; 
muziek[l , somij] := toonladder [a] ; 

end; 
end; 

{tweede rij opvullen} 
if k > 1 then 
begin 

for j := 1 to (k-1) do 
begin 
somij := i+j-1; 
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muziek[2,somij] := 1 
end; 
end; 

somik := i+k-1; 
muziek[2,soinik] := 0; 

i := i + k; 
if k = then i := aantaltijden+1 ; 



end; 



{laatste element tweede rij moet nul zijn} 
muziek [2, aantaltijden] := 0; 

{wegschrijven willekeurig muziekstuk} 

schrijfabc (muziek, bestandsnaammuzl) ; 

writeln(logf ile , '*******************************') 

writelndogf ile, '* Logfile TabuMusic *') 

writeln(logf ile , '*******************************') 

writelndogf ile) ; 

writelndogf ile) ; 

writelndogf ile, 'De willekeurig gegenereerde muziek is 

opgeslagen als ' , bestandsnaammuzl) ; 

writelndogf ile) ; 

writelndogf ile, 'Grondtoon: ', #9, grondtoon) ; 

writelndogf ile, 'Aantal maten: ' , #9, aantalmaten) ; 

writelndogf ile) ; 



{berekenen score willekeurig muziekstuk} 

show := 0; 

punten := doelfunctie (muziek) ; 

{aspiratieniveau initialiseren} 

for i :=1 to 2 do 
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begin 

for j := 1 to aantaltijden do 

begin 

aspiratie [i , j] := muziek[i,j] 

end; 
end; 

aspiratiescore := punten; 



writelnC '*******************') 

writelnC* TabuMusic *') 

writelnC '*******************') 

writeln; 

writelndogf ile, 'De score van het willekeurig muziekstuk 



bedraagt : ', punt en: 0:6); 

writelndogf ile) ; 

show := 1; 

punten := doelf unctie(muziek) ; 

show := 0; 

writelndogf ile) ; 

writelndogf ile, ' ') ; 

writelndogf ile) ; 

writelndogf ile, 'Er werd geoptimaliseerd met steeds ', iters, 

' iteraties : ' ) ; 

writelndogf ile, 'Lengte tabulijst melodie: ', #9,aantaltabu) ; 
writelndogf ile, 'Lengte tabulijst ritme : ', #9, aantaltabu2) ; 
writelndogf ile) ; 

writelndogf ile, ' ') ; 

writelndogf ile) ; 

writelnClogfile, 'WEGINGSCOEFFICIENTEN: >); 

writelndogf ile) ; 

writelnClogfile, 'Wegingscoeff icient Pitch variety: ', #9, 



wegingl :0:3) ; 

writelnClogfile, 'Wegingscoeff icient Dissonant intervals:', #9, 
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weging2:0:3) ; 

writelndogf ile, 'Wegingscoeff icient Contour direction: ', #9, 

weging3:0:3) ; 

writelndogf ile, 'Wegingscoeff icient Contour stability: ', #9, 

weging4:0 :3) ; 

writelndogf ile, 'Wegingscoeff icient Rhythmic variety: ', #9, 

wegingS :0 :3) ; 

writelndogf ile, 'Wegingscoeff icient Rhythmic range: ', #9, 

weging6:0:3) ; 
writelndogf ile) ; 

writelndogf ile, ' ') ; 

writelndogf ile) ; 

writelnClogfile, 'OPTIMALISATIE: '); 

writelndogf ile) ; 



writelnC 'Begins core : ' ,punten:0: 6) ; 
writeC'Bezig met verwerken. . . ' ) ; 

s := 1; 

while (s = 1) and (error) do 



begin 
show := 0; 

readlndnf ile, kies) ; 

if kies = 1 then 
begin 
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writelnClogfile, 'MELODIE: '); 

writelndogf ile) ; 

{STARTTIJD} 

start := Time; 

{tabu search melodie} 

show := 0; 

for d := 1 to iters do 

begin 

{tabu search: verbeteren toonhoogte noten} 

muziek := noottransponeren(muziek, punten) ; 

if d = iters then show := 0; 

punten := doelfunctie (muziek) ; 

show := 0; 
end; 

write (' . ') ; 
einde := Time; 

verschil := einde - start; 

DeCodeTime(verschil , uur, min, sec, ms) ; 
totaal := totaal + verschil; 

writelnClogfile, 'Score: ', punten:0 : 10) ; 

writelnClogfile, 'Beste score: ', aspiratiescore :0 : 10) ; 

writelnClogfile, 'Rekentijd: ', uur, ' uur, ', min, ' min, ', 

sec, ' sec, ' , ms, ' ms') ; 
writelnClogfile) ; 
writelnClogfile) ; 

end; 



if kies = 2 then 

begin 

{tabu search ritme} 

writelnClogfile, 'RITME: '); 

writelnClogfile) ; 

show := 0; 
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hulp := punt en; 

start := Time; 

for d := 1 to iters do 

begin 

muziek := ritmeverbeteren(muziek, punt en) ; 

if d = iters then show := 0; 

punten := doelfunctie (muziek) ; 

show := 0; 

end; 

write (' . ') ; 
einde := Time; 
verschil := einde - start; 
DeCodeTime(verschil , uur, min, sec, ms) ; 
totaal := totaal + verschil; 
writelndogf ile, 'Score: ', punten:0 : 10) ; 

writelndogf ile, 'Beste score: ', aspiratiescore :0: 10) ; 
writelndogf ile, 'Rekentijd: ', uur, ' uur, ', min, ' min, ', 

sec, ' sec, ' , ms, ' ms') ; 
writelndogf ile) ; 
writelndogf ile) ; 
end; 



if (kies <> 1) and (kies <> 2) then 

begin 

s := 0; 

end; 



end; 



writeln; 

writeln; 

writelnC'De beste score is: ', aspiratiescore : : 6) ; 

writelnC'De output is geschreven naar ' ,bestandsnaamlog) ; 

writelndogf ile, ' ') ; 
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writelndogf ile, ' ') ; 

writelndogf ile) ; 

writelndogf ile, 'Het GEOPTIMALISEERDE muziekfragment is opgeslagen 

onder ' , bestandsnaaminuz2) ; 

DeCodeTime(totaal, uur, min, sec, ms) ; 
writelndogf ile, 'Totale rekentijd: ', uur, ' uur, ', min, ' min, ' 

sec, ' sec, ' , ms, ' ms') ; 

writelndogf ile) ; 

writelndogf ile , ' EINDSCORE : ' ) ; 

writelndogf ile) ; 

show := 1; 

punten := doelf unctie(muziek) ; 

writelndogf ile) ; 

writelndogf ile, 'EINDTOTAAL: ' , #9, aspiratiescore : : 6) ; 

schrijf abc(aspiratie, bestandsnaammuz2) ; 

closednf ile) ; 

closedogf ile) ; 



end. 



E 



Voorbeeld van door Tabu- 
Music gegenereerde muziek 
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De score bedraagt 35. 



